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SYNOPSIS 

The topic of the study has been chosen from several 
considerations. In hard rock areas in India there, is a need 
for systematic investigations for groundwater development. 
Several of the drought-prone areas are situated in limestone 
terrain which if properly studied offers a potential source 
of groundwater. Limestone terrains contain typical features 
such as fractures, collapsed zones and solution cavities. 
Remote sensing offers a good scope for the identification and 
delineation of such features. Although studies on remote 
sensing for groundwater exploration have been carried out in 
India during last 10 years, most of them have been based on 
the visual interpretation of the imageries and the false 
colour composites that are directly obtained from National 
Remote Sensing Agency. There is a need for image processing 
of digital data and selective techniques for 


further 



enhancement of features that are of interest in groundwater 
exploration. To date, no such attempt has been made on a 
systematic basis for the carbonate terrains in India. 

Thus in the present work, both the aspects are 
dealt with. '■ Methodology and the sequence of operation for 
digital image processing are presented.* ^Enhancement 
techniques such as band ratioing, histogram equalization, and i 
edge enhancement as developed Iby the author and software for 
all the steps involved are provided. As a case study for the j 
application of this methodology, central part of the 
Chhattisgarh basin in Madhya Pradesh has been chosen. This | 
region supplies limestone to the Bhilai steel plant and many ■ 
other industries. Uith increasing mining and industrial I 

activity, demand for water has stepped up considerably in ' 

1 

recent years. Groundwater is the major source of supply. 
Regional field mapping has been undertaken by governmental 
agencies during the last six years. Uork on digital image 
processing and interpretation techniques contained in the 
study on the basin are reported for the first time. The 
methodology developed here is applicable to other drought- : 
prone limestone regions such as the formations of Cuddapah ; 
and Vindhyan supergroups in India. ; 

The study area covers around 2500 sq.kms. f The 
geological and drainage maps have been prepared on the basis 
of field work, aerial photographs and processed imageries. ^ In 
addition, measurements of static water levels for the pre- 



monsoon and post-monsoon periods have been made in larfee 
number of wells within the period pertaining to the analyzed 
satellite data, on the basis of which water level contour 
maps have been prepared. In the terrain delineated for good 
aquifer zones on the basis of digital image processing in the 
present work, pumptests were conducted on few selected wells 
and aquifer parameters were estimated. A few pumptests were 
also conducted on wells located outside the delineated zones 
for comparison with the productive aquifers. For making the 
groundwater picture complete, chemical analyses for limited 
number of samples of groundwater have been carried out. 

The work contained in the thesis is organized into 6 
chapters. In chapter 1, the topic has been introduced and the 
details of the study area are given. The preview of the 
Chhattisgarh basin with the significant contributions of all 
earlier workers is presented. The geology and characteristics 
of limestone and shale members within the study area are 
given in detail. In chapter 2 mathematical background of 
image processing techniques along with the methodology are 
dealt in detail. Several steps are involved in the pre- 
processing, image enhancement and image classification 
stages. The pre-processing stage contains all geometric 
corrections and data transfer for the study area. Conical 

present study 


orthomorphic projection has been used in the 
for geometric corrections as it is appropriate 


for a country 



like India with prominent north-south elongation. Utilizing 
the line numbers and pixel numbers generated, the study area : 
has been windowed from the tape an<5 processed for further- 
image processing work. The image enhancement includes several 
enhancement techniques in a sequence. Band ratioing forms the 
first step to achieve tonal contrast in the valley regions. 
Stretching forms the next step involving both the linear and 
selective stretching. The histograms initially plotted over 
the range of reflectance values ( 1.0 to 2 . 5 ) have been 
subsequently stretched to the entire medium which has 

resulted in considerable improvement of the resolution of the 
area. The lithological differences in limestone and shale 
members have been clearly brought out. Pseudocolouring 
technique was used to improve the morphological 
characteristics of three limestone members. Histogram 

equalization developed by the author has been used in | 

combination with stretching to enhance vegetational areas i 
that indicate the proximity of groundwater. A 5 by 5 filter | 
as developed in the present work has been used for edge 
enhancement to clearly bring out the edges of features like 
sinkholes and collapsed zones within the limestone members. ! 
The entire area has been classified on the basis of training 
sets with the help of Bayes’ classifier. All the i 

computational steps involved and the software developed for 

^ ' ' ■ ' f' 

all these techniques are presented in this chapter. | 

The lithological identification of the formations 



is given in chapter 3. On the basis of image enhancement 
techniques the area has been classified into three 
lithological units namely limestone, shale and alluvial 
tract. Uithin the limestone tract, three members have been 
demarcated based on the reflectance values, tonal 
characteristics and the soluble nature of the members. The 
details of the characteristics of all these formations 
relevant to groundwater are provided. 

In the fourth chapter, various morphological 
features, structural details and drainage patterns in the 
study area along with the land use classification are 
presented. The study area was classified into three distinct 
morphological units- upland tract, valley tract and moisture- 
stress tract. Areas with fractures and solutions features 
have been delineated using the techniques presented in the 
previous chapter. The drainage behaviour in the cavernous 
zones has been worked out. Classification of the study area 
into land use and land cover types has been carried out and a 
line printer map has been prepared. The accuracy involved has 
been justified on the basis of statistical analyses. The 
groundwater situation in general and the evaluation of the 
delineated aquifer zones in particular form the theme for 
chapter 5. The groundwater level contour maps for both pre- 
monsoon and post-monsoon periods are prepared and the areas 
of recharge as well as the possible directions of the 



T 

£roundwater flow in the study area are identified. The 

cavernous zone in the central part of the area acts as a sink 

for groundwater flow. Of the productive aquifer zones 

identified within the three liinestone members namely the 

Nipania, Pendri and Newari members, the Pendri member is 

identified as the most productive one with the t ransmissivity 

2 

(T) at places ranging to 1150 m /day. The productive zones in 

2 

the Newari member registered a maximum T value of 700 m /day 

2 

while the Nipania limestone has a range of 100 to 124 m /day 

for the same. Estimates of T from pumptest data conducted on 

wells outside the productive zones in all the three limestone 

members have revealed the T values to be of the order of 50 
2 

to 60 m /day. The details of the aquifer evaluation of the 
productive aquifer zones in the study area delineated on the 
basis of image processing have been presented. The variation 
of the groundwater quality as related to groundwater flow 
direction is also discussed in this chapter. The summary of 
the findings in terms of a comprehensive hydrogeological 
picture and the applicability of remote sensing for 
groundwater are presented in chapter 6. 



CHAPTER 1 


INTRODUCTION 

1.1 GENERAL 

The central region of the Chhattisgarh basin, 
chosen for the present work, has been a center of mining 
activity for the last three decades as it constitutes an 
important source of supply of dolomitic limestone to the 
Bhilai steel plant and also limestone for few cement 
factories. Groundwater is the principal source of water 
supply. Demand for water has been on the increasing trend 
with the increasing industrial activity in this region. 
Although sporadic attempts have been made to develop 
groundwater at certain chosen locations, no systematic effort 
has been made to locate favourable zones within this 
carbonate terrain. The terrain is flat with gently undulating 
topography. Areas with considerable groundwater storage exist 
in the region due to the recharge from surface drainage and 
the formation of solution cavities through selective 
dissolution of limestone. 

In the present work, delineation of such potential 
aquifer zones has been carried out through digital image 
processing techniques. To date, digital interpretation of 
imageries has not been attempted for formations in this 
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refiion by any of the earlier workers. The author, for the 
first time, has adopted an approach to the problem of the 
location of aquifer zones within the carbonate terrain 
through filtering and edge enhancement techniques. The 
software used in this work at different stages of operation 
has been developed during the present study. Systematic 
pumptests conducted in the delineated zones have confirmed 
the presence of such potential aquifer zones. 

1.2 CHHATTISGARH BASIN -A PREVIEU 
1.2.1 General 

The Chhattisgarh basin, one of the major sedimentary 
basins in India, is a crescent shaped intra-cratonic basin. 
It occupies an area about 33000 sq.kms. between the latitudes 
20*30’-22*20’N and longitudes 80 M 6 ' -83 * 46 ’ E , within Durg , 
Raipur and Bilaspur districts of Madhya Pradesh and adjoining 
areas in Orissa (Figs. 1.1 and 1.2) and has an annual rainfall 
of 110-132 cms,. The climate is tropical with temperatures 
fluctuating between 11*C and 47.5”C. 

The sediments mainly comprise of arenaceous and 
calcareous sequence attaining a total thickness of 550-600 m 
and unconf ormably overlying the Archean basement of granites, 
gneisses and metasediments. The formations exhibit varying 
strikes with gentle centripetal dips, except along the 
western and northern fringes where the basin has faulted 
margins against metavolcanics and granites. Part of the 
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southern boundary is also £neissic. Outliers of Gondwana 
formations represented by sandstones with basal conglomerate, 
occur in the eastern part. These formations are intensely 
folded and faulted. Disturbances along the northern margin 
are also evidenced from the occurrence of faulted ridge of 
Chandarpur sandstone near Karsia. In general, the shales and 
limestones exhibit flat topography with intermittent 
sandstone ridges. The southern and southwestern margins of 
the basin do not exhibit any signs of disturbances. 

The earliest reports on the geology of different 
areas in the basin date back to over a century (Meddlicot, 
1867; Blanford, 1870; King, 1885; Ball; 1897, Fermor; 1936). 
These investigations in general consisted of descriptions of 
field features and lithology. Due to the paucity of detailed 
information, the age correlations with other- stratigraphic 
sequences as reported by these authors were not accurate. 
Pascoe (1950) subdivided the sediments in the basin into two 
broad categories, the quartzo-conglomeratic Chandarpur series 
and overlying calcareous, clayey and marly Raipur series, of 
essentially a shale-limestone sequence, with the Kurnool 
group and Vindhyan supergroup of formations. 

Dutt (1964) on the basis of detailed field mapping 
in southern part of the basin has reported the first 
comprehensive stratigraphic succession of Chhattisgarh series 
(Table 1.1) and also indicated the possibility of correlating 
the same with sediments of adjoining Kurnool and Bastar 



TABLE 1.1 

STRATIGRAPHIC SUCCESSION OF CHHATTISGARH BASIN 
(After Dutt, 1964) 
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groups. In the study area located in the central part of the 
Chhattisgarh basin, the formations belong to the Raipur 
series. 

1.2.2 Raipur Series 

Schnitzer (1974) has subdivided the Raipur series on 
the basis of detailed investigations. He has reported the 
cyclic nature of the sedimentation within the same. The 
series has been divided into five cycles (Table 1.2). The 
formations within each cycle have been demarcated on the 
basis of occurrence of ripple marks and stromatolites as also 
the Ca and Mg contents of limestones. The formations of the 
first cycle are well exposed between Sarangarh and 
Seorinarayan while the second cycle formations occur to the 
west of this area. The Akaltara-Li lagar region is the type 
area for third cycle of formation. The area between Raipur 
and Bilaspur has fourth cycle of formations with the 
occurrence of fifth cycle of formations in the locations S 
and SU of Bilaspur. 

Based on the oligomictic conglomerate, coarse 
grained nature of the sandstone and ubiquitous cross-bedding 
(Fig. 1.3) the Chandarpur series underlying the Raipur series 
have been relegated to beach environment of deposition 
(Hurt i , 1 987 ) . The dominance of limestone and shales in the 
sequence of Raipur series is an indication of increase the 
depth of the basin during its deposition. The even lamination 
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TABLE 1-2 CYCLES OF SEDIMENTATION IN CHHATTISGARH 
BASIN (After Schnitzer ,1976) 

























































of bedded chert suggests its deposition in quiet water 
conditions. The existence of fine laminations in the 
limestone points to the lack of strong currents which is to 
be expected in a subtidal environment or sheltered lagoon 
(Hurti , 1987 ) . According to Murti, the micritic nature, 
thinness of laminations and general absence of sedimentary 
structures show that the limestone was formed under quiet 
water and low-energy conditions below wave base. The overall 
paucity of elastics is indicative of the peneplain nature for 
the provenance. Schnitzer (1974) inferred that the formations 
of Raipur series have been deposited in sub-tidal to tidal 
environments as each of the cycles starts with limestone and 
ends with shale. Minor uplift of the landmass by the 
epirogenic forces possibly released fine mud which ultimately 
resulted in the formation of shale. Dominance of carbonates 
in the sequence reflects prevalance of warm temperatures, 
absence of clastic supply, chemical weathering in the 
provenance and removal of CO 2 from sea water. Extensive 
occurrence of dolomites and dolomitic limestones indicates 
adequate availability of Mg in sea water. The pink shale of 
Gunderdehi formation shows that it was formed in an oxidizing 
environment. Dominance of the pink shale (Gunderdehi 
formation) reflects the humid and warm climate (oxidizing 
environment) that has prevailed in the region. The frequent 
interbedding of shale with the Pendri and Nipania limestones 
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FIG'1-4 REGIONAL CORRELATION ON THE BASIS OF CYCLIC SEDIMENTATION 





is of interest. The former is interbedded with olive green 
shale which is indicative of the deposition in a subtidal or 
oxygen-poor environment. Olive green and purple shale 
interbedded with purple limestone of Nipania reveals the 
nonuniformity of the conditions of the deposition. Purple 
shale intercalations reflect an oxidizing environment. 
Further, the basin might have fluctuated between oxidizing 
and reducing environments in the Nipania times. 

The regional correlation of the formations in the 
northern Chhattisgarh basin with the corresponding ones in 
Bastar, Jagadalpur and Bhima groups as also with the Kurnool 
group of the Cuddapah supergroup is presented in Fig. 1.4. 
Murti (1987) considered the formations in the Chhattisgarh 
basin to be equivalent to the formations of Rewa and Bhander 
series of Vindhyan supergroup (Table 1.3). Murti has given a 
comprehensive correlation of the stratigraphy of the entire 
chhattisgarh basin (Fig. 1.5). The dominance in the eastern 
and western margins of the basin with shale-limestone 
sequence in the central and southern regions is indicative of 
shallow nature of the basin towards the margins. The 
dominance, of sandstone in the eastern and western margins of 
the basin with shale-limestone sequence in the central and 
southern regions is indicative of shallow nature of the basin 
towards the margins. 

Estimation of the age of the formations in the 
Chhattisgarh basin has been carried out by several workers. 



STRATIGRAPHIC CORRELATION OF THE KURNOOL ,CHHATTISGARH 
GROUPS AND VINblHYAN SUPERGROUP 

Kurnool group Chhatt isgarh group Vindhyan supergroup 
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FIG.1-5 STRATIGRAPHIC CORRELATION IN THE CHHATTISGARH BASIN 




Uhi]e Dutt (1964) reported the sediment to be younger than 
1400 m.y, Schnitzer (1971) estimated the age to be between 
1000-600 m.y. Sarkar (1972) conceived the Chhattisgarh 
sediments to be younger than 600 m.y. Krenzer and co-workers 
(1977) assigned 700-750 m.y. for Chandarpur sandstone on the 
basis of K-Ar dating of the authigenic glauconite. Verma and 
co-workers (1977) on the basis of pal eomagnetic studies have 
concluded that the pole position for pink shale of Gunderdehi 
is intermediate between Rewa and Cuddapah sandstones. Since 
the age of Rewa formations is reported to be 800-900 m.y 
(Crawford , 1 960 ) and that of Cuddapah is 1220-1600 m.y 
(Crawford and Compston , 1 9 73 ) , the age of Gunderdehi formation 
on the basis of pole position is likely to be close to 1250- 
1300 m.y. 

1.2.3 Drainage 

The drainage network in the basin is a part of the 
Kahanadi system. The Mahanadi river has its course from SU to 
NE and drains through northeastern part of the basin 
(Fig. 1.2). Several of its tributaries, ephemeral and 
perennial in nature, drain across the basin. The prominent 
one is the Seonath river which is a result of the confluence 
of Tandula and Kharun in the central part of the basin. The 
Jamunia river joins Seonath at Rajpur. Contribution of 
surface drainage to groundwater is significant in limestone 
terrain in the study area. 


1.3 GEOLOGICAL SETUP OF THE STUDY AREA 


1.3.1 General 

Present work covers an area of around 2500 sq.kms. 
within the Raipur, Baloda Bazar, and Bilaspur districts in 
the central part of the basin. The formations in the study 
area (Fig. 1.6) pertain to Raipur series. Only members 
belonging to the Chandi and Tarenga groups are present. 
These formations correspond to the 3rd and 4th cycles in the 
sequence (Schnitzer, 1969). The stratigraphic succession and 
the geological map for the same are presented the Table 1.4 
and Fig. 1.6 respectively. 

The Raipur group of sediments overlies the 
Chandarpur group with an erosional disconf ormi ty . The 
disconf ormabl e nature of the contact is very clear between 
Birkoni and Tumgaon, where limestone has been deposited on 
the eroded surface of the sandstone. 

1.3.2 Chandi Formation 

The Gunderdehi formation is succeeded conformably by 
Chandi formation (comprising limestone members) with 
gradational contact. The formation is characterized by 
prolific growth of stromatolites (Fig. 1.7) occurring as 
extensive reef complexes in the area. Uithin the Chandi 
formation, sandstone pockets of localized occurrence have 
been reported by Dutt(1964). This sandstone, designated as 
Raipur sandstone at the base of Akaltara limestone sequence. 
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TABLE 1*4 STRATIGRAPHIC SUCCESSION OF THE STUDY 
AREA 






exhibits cross beddinfi with reversal of current directions 
(Fig. 1.8), indicative of an agitated environment. The Chandi 
formation can be broadly divided into three members 
(1) Nevari (2) Pendri and (3) Nipania on the basis of 
lithological considerations. 

1.3. 2.1 Newari Member 

The pink and purple coloured medium-grained limestone 
representing the Newari member conformably overlies the 
Gunderdehl formation with a gradational contact. In the 
vicinity of its contact with shale, the limestone is friable. 
Dolomitization is noticed along the joints within the same. 
Except in its proximity to the Gunderdehi formation, this 
limestone member is very massive and compact with indistinct 
bedding. Stromatolites and stylolites are quite common. The 
stromatolite columns are separated by calcareous mud and 
stand out prominently on weathering. 

1.3. 2. 2 Pendri Member 

The limestone member is gray in colour, fine-grained 
and is st romatol it ic . The limestone is characterized by 
frequent lithological variation to a greenish flaggy type and 
often with intercalations of shale. There is preferential 
dolomitization in this algal limestone. Stromatolites tend to 
be stunted, show branching and are ill-preserved because of 
the extremely fine lamination. The Pendri member corresponds 
to Bhatapara limestone of Schnl tzer C 1 96 9 ) . 
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1.3. 2. 3 Nipania Member 

Coarse to medium grained, this limestone member is 
purple with mottled appearance on the surface. Stromatolites 
with variety of species are abundant in the limestone. This 
member corresponds to Nandini limestone of Schnitzer (1969). 
The nature of stromatolites reveals the depositional 
conditions. The random orientation of the individual colonies 
of stromatolites is indicative of an agitational environment. 
Presence of coarser laminae in the stromatolites (Fig. 1.9) 
further confirms such an inference. The gradation of the 
stromatolites to finer laminae in the sequence is indicative 
of transition to a relatively quiet depositional environment. 

The members of the Chandi formation exhibit typical 
karst topography in the field. Differential weathering and 
structure has controlled the development of the landforms in 
the area. The collapsed zones (Fig. 1.10) on the surface in 
the field reveal the presence of sinkholes underneath. Karst 
morphology (Fig.l.(l) has a typical appearance of its own with 
a very thin and remarkably smooth patina of siliceous 
material on the karstic surface resulting in a relict fossil 
appearance on the karstic surface* 

Karsti f icat ion is structurally controlled. On the 
surface of an extensively exposed bed of limestone as in 
Kesili area, initially the swallow holes (Figs. 1.12 and 1.13) 
that have resulted from solution phenomena have rapidly 
coalesced giving rise to small channels which expand deeper 
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Channelised dints in Nipania limestone 









to form dints (Fig. 1.14). In the Raipur stroniatolitic 
limestone (Fig. 1.9), karst i f icat ion is developed along the 
contacts between the individual stromatolites as also along 
the contact of stromatolite with cement material, justifying 
thereby more prominent vertical karstif ication. Solution 
cavities (Fig. 1.15) are also noticed along the banks of 
rivers Jamunia and Seonath. Typical karst landforms are 
noticed in the study area in different forms, such as the 
sinkholes, lapies (Figs. 1.16 and 1.17) and caverns 
(Fig. 1.18). The sinkholes occur in different sizes with the 
depression varying between 0.5m and 8m . Sinkholes noticed in 
the study area are formed due to down-hole reaction. 
Uidespread outcrops of limestone occur having lapies with 
characteristic rugged, grooved, pitted, faulted and etched 
surface as around the Nawapara and Kesili villages. Cavities 
of different shapes and sizes are common along the banks of 
Jamunia (Fig. 1.15) and these act as conduits for groundwater 
movement . 

1.3.3 Tarenga Formation 

This formation is represented by a pink shale with 
interbedded chert. In the field, the occurrence of 
intraf ormational conglomerate is observed around Tarenga 
village. This formation corresponds to Patharia Umaria series 


of Schni tzer ( 1 9 6 9 ) . 





Karst landforia-Lapies in Newari limestone 




1.3.4 Hydrofieolo&y 

The karstic limestones in the present area 

constitute potential aquifers. Althoufih the limestone by 
itself possesses negligible primary porosity, with hardly any 
permeability when fresh and thus making it to be very poor 
aquifer. Karst i f icat ion due to potential dissolution with 
water has rendered the same into permeable complex having 
interconnected simple or complex conduit systems which 

facilitate free movement of groundwater in both confined and 
unconfined conditions. On the basis of field investigations 
in the shallow unconfined groundwater body of Raipur 
limestones, the open wells of moderate diameter with well 
developed cavities and interconnected conduit flow are 

capable of yielding 500-1000 m 3 /day with drawdown of 5 to 8 m. 
In such zones the individual wells have a specific 
capacity of 100 to 300 Ipm/m of drawdown and the cavernous 
aquifers with specific yield of 5 to 15 percent and hydraulic 
conductivity of 5 to 14 m/d. 

Published work on the hydrogeological and the 

aquifer characteristics of the Raipur limestone is scant in 

literature. Adyalkar and co-workers (1975) have reported the 

2 

values of transmissivity (T) to be between 40 m /day and 
2 

600 m while permeability variations are in the range 
of 0.9 m/day to 5 m/day. Chowdhary (1973) has discussed in| 
general the ground water provenances in Madhya Pradesh. 

no published reports are available on the aquifer 

2 ^ ! 


However , 



characteristics and the well performance in the different 
formations in the Chhattisgarh basin. 

1.4 SCOPE OF THE PRESENT STUDY AND ITS ORGANIZATION 

Digital image processing of satellite data plays an 
important role in the groundwater exploration and development 
in carbonate terrains. In India, as evident from scanty 
published literature, no attempt has been made for detailed 
study in this direction. The work reported herein deals with 
the development of the methodology and related software for 
identification of collapsed and cavernous zones and also 
sinkhole-prone areas that are potential zones of groundwater 
storage within carbonate terrain. The region from the central 
part of the Chhattisgarh basin is chosen as a case study for 
investigations. The demand for groundwater is ever increasing 
in this area with increase in mining and industrial activity. 
There is a need for systematic location of potential aquifer 
zones in the limestone terrain in this region and no work has 
been undertaken so far in this direction. Thus, the present 
study, an attempt for the first time in India, is undertaken 
to evolve suitable methodology and software for digital image 
processing for use in location of aquifer zones in carbonate 
terrain. Inferences arrived at, on the basis of digital image 
interpretation, in the study area have been confirmed by 
extensive ground checks through detailed field work. 

The material contained in the thesis is organized 

a9 



and presented under several chapters. In chapter 1, the topic 
is introduced with review of existing information pertaining 
to the Chhattisgarh basin in general and the study area in 
particular. The scope of the present study also outlined. In 

chapter 2, the methodology adopted in the present study is 

detailed with all the relevant mathematical background. Edge 
enhancement techniques developed in the present work are 

presented in detail. In chapter 3, the lithological 

identification within the limestone members and shale through 
enhancement techniques have been described. The morphology, 
structural features of limestone and their role in the 
formation of the caverns and collapsed zones together with 
the classification results of various landuse are given in 
chapter 4. Chapter 5 deals with the groundwater situation in 
general and the evaluation of the cavernous limestone 

aquifers in particular. In this chapter, analyses of the 
time-drawdown and time-recovery data for the pumptests 

conducted in the productive areas delineated on the basis of 
image processing interpretation are presented. The aquifer 
parameters of the limestone aquifers in the study area are 
estimated. Summary and conclusions are presented in chapter 6. 



CHAPTER 2 


IMAGE PROCESSING TECHNIQUES FOR 
HYDROGEOLOGIC INVESTIGATIONS 


2.1 GENERAL 

Remote sensing is the use of reflected and emitted 
energy to measure the physical properties of distant objects 
and their surroundings. The sensors in the range 0 . 4p.m-l . l;im 
of the electromagnetic spectrum (Fig. 2.1) are considered 
here . 

Photographic sensors still remain the primary tool 
of hydrogeologists. Black and white aerial photography is 
used to determine some of the parameters such as boundaries 
of the rock types and drainage pattern. Thermal sensors (8- 
14juim) have been used for several hydrological studies in the 
past. Stephenson (1 9 74 ) and Coker (1970) are two of many who 
found thermal imagery as a valuable tool in geohydrological 
studies, particularly with respect to the g-roundwater 
movement . 

In the present work, mult ispectral data have been 
used to delineate sinkholes and fracture zones which act as 
avenues for groundwater in carbonate rocks. Recourse has also 
been taken to aerial photographs for checking the drainage 
and moisture-saturation zones. Robinove (1965, 1966, 1968) 
has periodically reported the progress of remote sensing in 
the field of hydrology and hydrogeology. His evaluation of 
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FIG- 2-1 ELECTROMAGNETIC SPECTRUM 



remote sensors as hydrologic tools is presented in Appendix-I 
Information derived from these techniques is 
valuable in the study of carbonate terrains. Moreover, 
because large areas can be examined in a very short period of 
time, remote sensing technology can be viewed as a time 
saving tool . 

LaMoreaux and Uilson (1980) have reported the following 
uses of mul t ispectral data connected with hydrological 
studies in carbonate terrains: 

1. Inventory of sinkholes. 

2. Monitoring the sinkhole development. 

3. Mapping sinkhole alignments. 

4 . Invest igat ing the relationship among the sinkhole 
development, fracture traces, lineaments, and groundwater. 

5 . Del ineat ing collapsed zones. 

6. Mapping the surface karst springs. 

7. Locating potential water well sites. 

8. Mapping the regional geological structures. 

Moore (1980) has presented the various steps 
involved in the aquifer mapping using remote sensors. 
Details of the same are provided in the Appendix-I I. 

In India, studies on remote sensing have been 
reported since last 10 years covering agriculture, forestry, 
land use and mineral surveys. However, the work pertaining to 
groundwater investigations was limited and largely based on 
the visual interpretation of imageries. Santosh Kumar and 
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Ti mven^adachar i (1981) have presented an integrated approach 
for groundwater studies in parts of central and southern 
India on the basis of geoinorphology and hydrology information 
in relation to groundwater occurrence and movement. 
Chaturvedi and co—workers (1983) have evaluated the 
groundwater potential for the terrain in the Bundelkhand 
region of Uttar Pradesh based on Landsat imageries. 
Tillaigovindara jan et.al (1983) have evaluated the 
hydrogeological conditions in the Hadurai , Tirunelvel li , 
Ramanathapuram and Kanya Kumari districts of Tamil Nadu. The 
hydrogeological and lineament maps were prepared on the basis 
of the Landsat imageries (bands 4,5,6, and 7) and false 
colour composites. Studies have also been reported on the 
application of Landsat imageries for groundwater exploration 
in the state of Gujarat (Baldev Sahai and co-workers , 1 983 ) . 
The hydrogeological maps for the semi-arid hard rock region 
have been presented using Landsat imageries. Balakrishnan 
(1986) reported in detail the different aspects of the water 
resources development and management and has presented a 
summary of work done through remote sensing. As already 

indicated, in all these cases visual interpretation of 

, ■ 2 ' . 

imageries has been reported. Working with Landsat and M S in 

the Tirupati area of Andhra Pradesh. for the groundwater 

exploration, Sharma and co-workers (1985) have concluded that 

2 

the band ratio images of Landsat and M S as well as density- 



sliced and thermal pictures have better interpretation 
capabilities over the other remote sensing formats. They have 
delineated level~I and level-II land use and land cover 
categories for the Tirupati region. Rampal (1985) has given 
procedures for estimation of water table depths using digital 
reflectance data. Data for Kanpur region was used for 
classification in terms of the groundwater table depth and 
land use. A line printer map was also presented. During the 
last 5 years in India, efforts are concentrated on the study 
of groundwater conditions in limestone terrains, since 
initial invesigat ions have indicated a good potential for the 
same in extensive carbonate terrains particularly in the 
drought “prone areas in various states. No work has been 
reported to date in India on the digital image processing for 
aquifer delineation in limestones. 

It is easily seen from the mat erial presented above 
that there is a wide scope for remote sensing applications in 
hydrogeology. In order to be able to utilize or even propose 
a proper remote sensing system as a hydrogeological tool, one 
has to not only have a complete understanding of the physics 
of the electromagnetic spectrum but also of its interaction 
with hydrologic features. The ultimate objective of this 
study is to explore the use of remote sensing system to 
identify the source areas. 

2.2 ACQUISITION OF REMOTE SENSING DATA 

Space remote sensing began in earnest during the 



REMOTE SENSING 
SATELLITES 



FlG.2-2 TAXONOMY OF REMOTE SENSING SATELLITES 


period 1946 to 1950. The taxonomy of remote sensing 
satellites is shown in Fig-2. On July 23 of 1972 , NASA 
launched Earth Resources Satellite (ERTS), which produced 
moderately high resolution imagery of the cloud-free Earth. 
Just prior to the launch of ERTS-B on January 25 in 1975, the 
ERTS program has been officially redesignated by NASA as the 
’’Landsat Program” . 

In the present work, data from Landsat-4 have been 
used. The Landsat-4 moving in a Sun-synchronous orbit at an 
altitude of 705 kms with 16 days repetitive coverage having a 
field of view of 14.93* carries a four channel mult ispectral 
scanner and Thematic Mapper. On the NASA Landsat-4 coverage 
index map, the study area is located by path number and row 
number as 142-045. The Landsat data corresponding to the 
aforesaid scene was obtained in the form of CCT (1600 BPI) 
containing reflectance values. 

2.2.1 Mult ispectral Scanner System (MSS) 

The mult ispectral scanner system records four images 
for a scene in green, red and two near-infrared wavelength 
bands (Table. 2.1) which are identified by the channels they 
occupied in the satellite telemetry system as 1,2,3 and 4. 
Each frame of the imagery covers an area of about 185 by 185 
sq.km, on a scale of 1:1000,000 and has a ground resolution 
of 79m by 79m. Due to the 10 percent overlap, the pixels of 
the MSS image are rectangular instead of square in shape 



table: 2-1 



Landsat MSS 

bands 


Landsat 

1,2 & 3 

Landsat 

4 a 5 

Bandwidth 
( jjm) 

Waveband 

name 

4 

1 

0-5- 0-6 

Green 

5 

2 

0-6- 0-7 

Red 

6 

3 

0 

1 

o 

CP 

Near Infrared 

7 

4 

o 

1 

CO 

• 

o 

Near Infrared 


S8 





One of SIX fines 
imaged simuUaneousty 
in each scan 


3 240 Pixels 
per scan line 


V / 

2340 Scan lines 


FIG.2-3 CONFIGURATION OF LANDSAT MSS 
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FIG.2-4 THE SAMPLING PATTERN OF THE LANDSAT 
MULTISPECTRAL SCANNING SYSTEM 
(ModHied from NASA 1976) 



(Slater, 1979), consequently the ground coverage of each 
pixel being 56m by 79m (Fig. 2. 3 and 2.4). 

A complete channel of Landsat MSS imagery contains 
2340 scanlines and 3240 pixels per line with a total of 
around 7.5 million pixels per scene of 4 channels and the 
entire data recorded was recorded on computer compatible 
tapes (CCT ) . 

2.2.2 Thematic Mapper (TM) 

The Landsat 4 and 5 carried Thematic Mapper, which 
has higher resolution than that of MSS. It records 256 
radiance levels in 7 wave bands with a spatial resolution of 
30m in 6 of them and 120m in the seventh one, in contrast to 
MSS that records 64 radiance levels in four wave bands with 
spectral resolution of 79m. The details of wave bands are 
indicated in Table 2.2. The Landsat-4 TM data was not 
available at the Indian Landsat Earth Stat ion( I LES ) , located 
60 Km from Hyderabad, due to the failure of the X band radio 
channel (in Feb 1983) that is responsible for the 
transmission of TM data directly to ground stations. However, 
the MSS data from Landsat-4 is not affected since a separate 
S-band radio channel is involved. 

Along with the satellite MSS digital data, the 
Survey of India toposheets of the study area on 1:250,000 and 
1:50,000 scale were also used. 
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The wave bands of Thematic Mapper 
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table: 2-3 

Coordinates of the Imagery 
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2.3 PRE-PROCESSING OF LANDSAT MSS DATA 


The digital data of Landsat-4 MSS obtained on a CCT 
was processed on fourth generation computer system to 
delineate and map the land use. The CCT was recorded on a 32- 
bit computer system, which is not compatible to the DEC-1090 
system (existing at the Computer Center of I.I.T, Kanpur) as 
the latter is a 36-bit configuration and hence necessitates 
conversion of the original CCT. Uhen the original tape is 
run on DEC-1090 system, it was reading 4 bytes and 4 bits of 
the next word, thus tampering the data. This problem was 
overcome by converting the original data tape to DEC 
compatible mode, keeping dummy bytes after each computer word 
and copied on another tape. For this work. OMEGA computer 
system of CAD Center of I.I.T, Kanpur has been used. 
Subsequently, the original tape data was used in other 
computer system, PDF 11/24 (which has an image processor) for 
the image processing work. The software involved for use on 
all these systems was developed by the author during the 
present study and has been used for processing the data. 
2.3.1 Band Separation 

The CCT (of 27th Dec 1984) used in the present study 
was in BIL format (Fig. 2. 5). The data for each image were 
band separated and kept in different files. This operation 
has been carried out on the OMEGA system using software in C— 
language as developed by the author. 



2.3.2 Uindowina 

The CCT of one Landsat MSS scene represents over 
34,000 sq.km of the Earth’s surface. In the first stage of 
operation, the area of interest was extracted, pertaining to 
bands 1,2,3 and 4, from the complete scene with the help of 
line and pixel number addresses as calculated from the 
knowledge of scene corners. The data on the test site 
covering 480 scan lines and 715 pixels have been windowed. 
The boundaries of the windowed data are as follows : 


Starting 

scan 

line 

--1036 

480 scan lines 

Ending 

scan 

line 

— 1516 


Starting 

pixel 


-- 644 

715 pixels 

Ending p 

ixel 


--1358 



The precise demarcation of the study area has been 
confirmed with the help of the Survey of India toposheets. 
The coordinates and the respective constants of the entire 
scene are presented in Table 2.3 
2.3.3 Histograms 

Histograms reveal the frequency of various digital 
numbers present in the data. Each digital number (DN) 
represents the reflectance of each pixel area, on a pre- 
assigned gray level scale. The gray level scale may be on 0- 
63 (six bit number), 0-127 (7 bit number) or 0-255 (8 bit 

number). The frequency distribution is of immense use in the 
image enhancement process as it aids in understanding the 



spectral contrast . 


2.3.4 Cartographic Projection for Mapping Landsat MSS Data 

The Landsat series of satellites map geodal surface 
of the Earth on to plane. These images therefore show various 
geometric distortions. All satellite motion-related errors 
have been corrected at the time of generating the Landsat 
imagery itself. However, errors due to non-linearity of 
ground track of satellite are not corrected and these produce 
images which cannot be compared directly with cartographic 
maps. In order to bring Landsat images on par with 
cartographic maps, it is necessary to transform the line 
number and pixel number of digital images into Cartesian 
coordinate system of latitude and longitude. Landsat orbit is 
inclined at an angle of approximately 81* to the equatorial 
plane of the Earth. The MSS aboard Landsat scene has a long 
strip of 185 km width which is oblique to meridional plane. 
It is therefore desirable to use a conformal map projection 
based on North-South strips. Conical orthomorphic projection 
has been adopted in the present work as it is appropriate for 
countries like India with prominent north-south elongation. 

2.3.5 Conversion of Geographical Coordinates to Conical 
Orthomorphic Coordinates 

This requires the availability of accurate maps of 
the image area and image identifiable ground control points 
such as highway intersections and small water bodies. 
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FIG.2-6 CORNER CO-ORDINATES OF THE IMAGE.-AND 
CONVERSION TO LINE NUMBER AND PIXEL 
NUMBER 


4 ? 




For conversion, the precise longitudes and latitudes 
of the four corners of the imagery (Table 2.3) are obtained 
from Survey of India toposheets on 1:250,000 scale. The 
procedure involved is as follows : 

The origin ’ 0 ’ ( Fig . 2 . 6 ) , chosen just outside the 
imagery edges, (Fig. 2.6) has geographical coordinates (Bq, 
Lq), where ’Bo’is the latitude and ’ Lq ’ is the longitude. Let 
(Bp, Lp) be geographical coordinates of any point on the 
imagery. The conical orthomorphic coordinates are 


^, = (P - m’)SinT (2.1) 

X = ai’ + X tan(T/2) (2.2) 

where .. 

T = L sin Bq 

where L = Lp - Lq 


i.e. the difference between the longitudes at any point ’p’ 
on the imagery and the origin ’0' 

3 4 5 2 

m m tanBg m (5+3tan Bg ) 

m’= m + + ^ + 5 

6Ro Ng 24 Rg Ng 120 Rg Ng 

.... (2.3) 

Rg = Radius of the curvature of the Earth spheriod at 
P (meridional) 

a (1 - e^) 

Rq = ^ (2.4) 

Tf(l -e sin Bg) 

a = major axis of the reference ellipsoid (6377277.6m for 
Everest ellipsoid) 

b = minor axis of the reference ellipsoid (6356075.0m, 

for India) 
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eccentricity 


2 
e 

e^ = (1 -b*/a*) = 1/300.8 (2.5) 

Nq = Normal to the surface at Bg (radius of curvature in 
the direction perpendicular to the meridian) 


No 


a 

2 . 2 „ . 1/2 

(1 - e sin Bq) 


( 2 . 6 ) 


Rjj, = mean radius of the Earth at mean latitude 


m = (Bp- Bo) (2.7) 

p = Nq cot Bq (2.8) 


Using the above equations, the orthomorphic coordinates (X,Y) 


are calculated. 

2.3.6 Transformations 

The conical orthomorphic coordinates (X.Y) can be 
transformed to the Landsat imagery coordinates (x,y) using 
the equations 

X = Ai X + Bi y + Cl (2.9) 

Y = A 2 X + B 2 y + C 2 (2.10) 

where 

X = Line number 
y = Pixel number 

and 


Ai , A 2 , A 3 , A 4 , A 5 , A 5 are constants of transformations. 

The values for the six constants can be determined 
from a minimum of 3 ground points which are the three 





corners of the imagery, whose coordinates (x,y) are computed 
from the earlier equations for their geographical 
coordinates. The values of X,y obtained are presented in 
Table 2M The line number and pixel number of these three 
points are also known since they are corners of the imagery. 
For evaluating the six unknown constants we need to solve six 


simultaneous equations CEqns.2.11 to 2.16). 

Xi = Ai XI + Bi yi + Cl (2.11) 

Yi = A2 X2 + B2 y2 + C2 (2.12) 

X 2 = Ai X 2 + Bi y 2 + Ci (2.13) 

^2 =-^2X2 + B2 y 2 ^2 ( 2 . 14 ) 

X 3 = Ai X 3 + Bi y 3 + Cl (2.15) 

^3 = ^ 2 x 3 + B 2 Y2 ^2 (2.16) 


These six simultaneous equatibns were solved using Gauss 
elimination method and the constants Ai , A 2 , Bi , B 2 , Ci , C 2 
were determined. 

The validity of the solutions was checked at 4th 


corner using the equations 

X4 = Ai X4 + Bi y4 + Cl ( 2 . 17 ) 

Y4 = A2 X4 + B2 y4 C2 ( 2 . 18 ) 


The line number and pixel number (Xp.yp) of any 
point on the imagery can be determined by solving equations 
2 . 1 9-2 . 20 using the values of six constants already obtained 
(Rampal , 19 7f5) . 



TABLE 2.4 


COORDINATES AND ORHTOMORPHIC COORDINATES FOR THE LANDSAT-4 
IHAGERY OF CHHATTISGARH AREA.M.P 


20 .718126 
22.325477 
22.655706 
21 .044899 


83 .066590 
83 .438342 
81 .732534 
81 .378365 


X 

X 

X 

X 


215142.585862787705 Y = -29983.33940945559 
251112 . 531892367146 Y = 149459.03199593853 
75274.108094912138 Y = 184484.91900469838 
39304.712218415477 Y = 5043.94876056306 






= ^1 Xp + Vp + Cl (2.19) 

Yp = A2 Xp + B2 Yp + C2 (2.20) 


The six constants of the imagery of Chhattisgarh 
Basin are presented in Table Z.t. 

2.4 COMPUTER PROCESSING OF THE LANDSAT-4 DATA 

Through computer processing of digital data , valuabl e 
information is obtained. The three properties of Landsat data 
which are important for computer processing are the 

mult ispectral nature of the sensors, the presentation of the 
data in the image format and multitemporal coverage of the 
orbits. Number of steps are involved in the computer 
processing of remotely sensed data. 

The methodology of the image processing techniques 

used in the present study is presented in the subsequent 

pages within this section. 

2.4.1 Image Restoration 

This is the first stage in any image processing 

sequence. . It was used to restore distorted image data to a 

more meaningful representation. It includes both radiometric 
correction and geometric restoration and correction. The 
radiometric corrections are done while preparing the computer 
compatible tape where as the geometric corrections are to be 
done to make it. conformal with base maps. 

2. 4. 1.1 Geometric Restoration and correction 

The aim of the geometric restoration and correction 
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TABLE 2.5 


CONSTANTS OF TRANSFORMATION FOR THE LANDSAT-4 IMAGERY OF 

CHHATTISGARH AREA.M.P 

A1 = - 15.37834371515483 

B1 = 54.28787394037448 

Cl = 75235.19856468876090 

A2 = - 76.71755874202162 

B2 = - 10.81379684548759 


C2 


184572 . 45036062989400 





is to make an image conform to a prearranged scheme. It can 
be done through resampling. Image resampling involves the 
reformation of an image onto a new base. This is achieved by 
using features that are common both in the image and the new 
base (Bernstein, 1976). These features which are termed 
ground control points (GCPs), are chosen to be in sharp 
contrast to their surroundings such as the edges of water 
bodies and airport runways (Bernstein, 1976; Benny, 1983). 
The ground control points are located on the image by their 
’x’ and 'y' coordinates and on the new base by their latitude 
and longitude. The functional relationship ( f and f 2 ) 
between image 'x' and 'y' and latitude and longitude are 
determined by least squares regression as follows : 


X = fi (latitude and longitude) .....(2.21) 

y = f2 (latitude and longitude) ......(2.22) 

where 


(x,y) = Image coordinates (column, row) 
and f^ , f2 = transformation relationships. 

Four stages are involved in the use of these two 
equations to geometrically transform an image : 

(1) A geometrically correct geographical grid is defined in 
terms of latitude and longitude. 

(2) The computer proceeds through each cell in this 
geographical grid and transforms the latitude and 
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lonfiitude values into values of x,y (Eqs.2.21 and 
2.22), which become new address in the image pixel. 

(3) The appropriate pixel value is transferred from the 

image data set to its new location in the geographical 
grid . 

C4) This process is repeated until the geographical grid 
is full at which point the image has been 

geometrically corrected usually to an accuracy of 
+ 1 pixel (Bernstein, 1 9 7 6 ) . T-yv 
There are three methods available for this 

resampling procedure. 

(1) Nearest neighbour method. 

(2) Bilenear interpolation. 

r 

(3) Cubic convolution. 

In the present study, nearest neighbour method is 
adopted as it is computationally simple and avoids altering 
pixel values. 

f(x) = fCxj^) for the X£ closest to x 
2.4.2 Image Enhancement 

Image enhancement aids in highlighting the features 
of interest by improving the visual appearance of the 
image, as also the interpretabil ity of an image by increasing 
apparent contrast between the features in the scene .Several 
steps involved in achieving the image enhancement are as 
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follows : 



0 

FIG-2-7 PRINCIPLE OF COI 
ENHANCEMENT 
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2. 4. 2.1 Contrast Stretching 

Contrast stretching is one of the most widely used 
enhancement techniques in image processing for lithological 
mapping . Sine e the Landsat MSS is designed to accommodate a 
wide range of scene illumination conditions , from poorly lit 
arctic regions to high reflectance desert regions, the pixel 
values in the majority of the scenes occupy a relatively 
small portion of the possible range of image values. If the 
pixel values are displayed in their original form, only a 
small range of gray values will be used , resulting in a low 
contrast display on which similar features might be 
indistinguishable. Contrast stretching is the process of 
redistributing the brightness levels in a picture to utilize 
the entire dynamic range of the display device (Fig. 2. 7). 
Several linear and non-linear algorithms are available which 
can be used for this technique . Depending upon the probability 
density function of the scene and identifying the areas of 
interest, a suitable algorithm can be chosen to maximize the 
contrast between those areas (Siegal and Gillespe, 1979). A 
linear stretch is generated by mapping each gray level 
encountered in the image, pixel by pixel, and then separating 
the pixel gray level as much as possible to fill the complete 
dynamic range of the display medium. The contrast stretch 
algorithm (Arl ington , 1979 ) utilized in the present work is as 
follows : 




Min 


X 265 



where 


G 


i j 

Max - 


Min 


(2.23) 


Gj-j = input pixel values line ’ i ’ and sample 'j* 

G'ij = corresponding output pixel value after 
the contrast stretch. 

(Min, Max) = The range of Min-Max brightness of a 

particular image 
2 ,4. 2. 2 Histogram Equalization 

A technique for obtaining uniform histogram is 
known as histogram equalization (Gonzalez and Uintz,1977). In 
the case of digital image processing, the gray levels are 
discrete values. Since original image is rescaled, this type 
of technique is useful to know in detail the darker regions 
of the image. Mathematically, it can be explained as follows: 

Let ’r’ represent the gray level of the pixels in 
the image to be enhanced. The pixel values have been 
normalized so that they lie in the range 
0 ^ r ^ 1 

with r = 0 representing black and r= 1 representing white in 
gray scale. For any 'r^ in the interval [0,1], attention is 
focused on transformations of the form 

S = TCr) ...(2.24) 

which produces a level 's’ for every pixel value 'r' in the 
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original image. It is assumed that the transformation 
function given in Eqn.2.24 satisfies the conditions 

(a) TCr) is a single value monot oni cal ly increasing in 
the interval 0 ^ r 1 , and 

(b) 0 f TCr) i 1 for 0 <: r ^ 1 

Uhile the first condition preserves the order from black to 
white in gray scale, the second one generates mapping that 
is consistent with allowed range of the pixel values. Keeping 
these conditions in mind, a transform function has been 
considered in the present study which is represented as 


s = TCr) = 


I p C^) 
0 


0 i r i 1 


. C2 . 25) 


where ’w’ is a duflimy variable of integration. The extreme 
right side of the equation . is recognized as ^ cumulative 
distribution f unction(CDF) of ’r’. The above concept has been 
formulated in discrete form for digital image processing. 

For gray levels which assume discrete values, we 
deal with probabilities given by the relation 

PpCr^ ) = 0 < r < 1 

n 

k = 0,1, 1-1 

...... (2.26) 


where L is the number of levels, Pr-Cri^.) is the probability of 
the k th gray level, n is the number of times this gray 
level appears in the image, and ’n’ is the total 
pixels in the image. 


number of 



The discrete form of Eqn.(2.25) is given as 


®k = TCr^) 



k 

= E pj. 
j = 0 


(rj) 


k 


0 i ry. i 1 

0,1 L - 1 


C2 . 27) 


Software for this has been developed by the author which has 
been used in the detailed study of densely vegetated collapsed 
zones. Some system routines (PDF 11/24) have been used to make 
it interactive . Results of the same are presented in the 
subsequent chapter. 

2. 4. 2. 3 Band Ratioing 

Band ratioing is useful in enhancing minor 
differences between geological formations. Ratio images 
minimize the first order effects due to variations in the 
brightness arising' from the slopes of the formations . Rat io 
images display the differences between reflectance spectra of 
bands in one image. Band ratioing is a processing technique 
in which two spatially registered spectral images are 
divided , pixel by pixel. The resultant data are rescaled to 
use the full dynamic range of display using contrast 
enhancement techniques. 

Mathematically, 

DN’ = X DNi (2.28) 

DN2 
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wher e 


DN’ is the ratio of digital number from channel 1 and 
channel 2 . 

The DN is bright when DN 2 > DN^ and dark if DN 2 <DNi.The 
ratio pictures have the effect of suppressing the details of 
features like the topographic effects or vegetational effects 
while emphasizing colour boundaries. This property has made 
ratio pictures quite useful in geological applications as 
they exaggerate subtle colour differences in a scene. 

2. 4. 2. 4 Edge Enhancement 

Edge enhancement filters are high-pass filters 
and their effect is to enhance or boost edges. This technique 
enables enhancement of subtle gray scale variations. Uhen the 
values of adjacent pixels vary by a more than predetermined 
threshold value, the interface is marked on the image display 
with a contour. It is achieved through spatial filtering. 
Spatial filtering is of two types. In one, the procedure is 
undertaken on the transformed image frequency data in what is 
known as "frequency domain” and the other one is known as the 


"spatial domain”. In 

the 

present 

v'ork, spatial 

domain 

filtering is adopted. 

The 

filtering 

in spatial 

domain 


involves passing a filter of N X N pixels over an image and 
converting the pixel of this so called "kernel” to some 
measurement of the total array. As this technique involves 
the rolling together of all DN within the pixel array, it is 
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Output at (4,3) 
multiply and sum 

= h(f32 + h2f33 + h3f34 
h4f42 + h5 f43 + h5 f44 
+ h7f52 + h8f53 + h9 f54 




Output at (4,4 ) 

= multiply and sum 

= h| f33 + h2f34 1^3 f35 

^4t43 ^5 W ^6^4 5 

*^ 7^53 " ^^ 8%4 


FIG.2*9 DIGITAL CONVOLUTION 



often known by a descriptive term as "Convolution ’ ( Kim and 
St r int z i s , 1 9 8 0 ) . The convolution is the movement (pixel by 
pixel and line by line) of a filter window through an image. 

Mathematical representation of the convolution is as 
follows : 

A system ’S’ is considered a black box with an input f(x) and 
output gCx) = S(fCx)). 

fCx) > S > gCx) (2.29) 

In this case , ’ f (x) ’ is the original image , represent ed in one 
dimension for simplicity , ’g(x) ’ is the filtered output image, 
and ’S’ the filtering operation . Then the output of the filter 
can be expressed as 


p 


fiCx) = 

where hCf) is the 
response. In digital 


fCt) h(x 

point 
case , the 


“ t ) dt 

spread function 
integral becomes 


..... (2 . 30) 

or impulse 
summation : 


+ O0 

gCi) = E f (k) hCi-k) .... (2.31) 

— oo 

and, although the limits of the summation are infinite, the 
function ’h’ is usually zero outside some range as indicated 
in Fig. 2.8. If the range over which ’h’ is non-zero is (-w to 
+w) , the equation (2.31) may be written as 

i+w j+w 

g(i,j) = E E f(k,l) h (i-k,j-l) ..(2.32) 

k=i-w l=j-v 

Thus the output g(i.j) at a point i,j is given by a weighted 
sum of input pixels surrounding i,j where the weights are 


g4 


given by h(k). To create the output at the next pixel i+1 the 
function h(k.) is shifted by one, and the weighted sum is 
recomputed. The full output is created by a series of shift- 
multiply-sum operations (termed as digital convolution) shown 
in Fig. 2. 9. The values of 'h’are referred to as the filter 
weights, the filter kernel, or filter mask. 

The high-pass filtering is of two types-fl) gradient 
filters and (2) laplacian f liters . These are defined for a 
continuous function fCx,y) of two variables as 
Gradient 


6f 


6 £ 


Vf 


Laplacian 

2 . 
^ f 


i + 
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(2 . 34) 


fix fiy 

where i and j are unit vectors in the x and y directions, 
(a) Gradient 


The gradient operator applied to a continuous 
function produces a vector at each point whose direction 
gives the direction of maximum change of the function at 
that point and whose magnitude gives the magnitude of this 
maximum change. Generally Sobel operators are used for 
Landsat data in this type. Few typical gradient operators are 
presented in Fig. 2. 10. The summation of the elements in each 
filter equals to 'O’. 
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(b)Laplacian filter : 


filter has been used in the present work to 
enhance the sinkhole and collapse prone areas in the study 
area. This filter is also computed by convolving the image. A 
sharpened edge without regard to edge direction is obtained 
through this process. 

riathemat i cal ly the filter can be derived as 

Continuous Digital 

fCx) vCi) 

f’Cx) v’Ci) = v(i-l) 

f^Cx) =V fCx) v”Ci) = vti) - vti-1) 

= [vCi)-v(i-l)]-[vCi-l)-vCi-2)] 

= v(i-2)-2v(i-l) tv(i) 

= Cl-2 1) (vCi-2) v(i-l) v(i)) 

C2.35) 

giving the convolution mask (1 ~2 l).In this form the 

laplacian at 'i' is computed from the values centered about 
’i~l’. To keep the laplacian symmetric, it is normally 
shifted and given at ’ i ^ as : 

(1 ~2 1) vCi~l) vCi) v(i+l) 

This is a one-dimensional digital laplacian filter. Based on 
the same principle, in the present work a two-dimensional 
filter has been used. Few examples of the two-dimensional 
laplacian filters have been presented in Fig. 2. 11. 

Though the digital image processing technique gives 
more details in recognizing the source areas, it is restricted 
to a spatial constraint of the processor window. Due to this 
reason a line printer map has also been presented to 
represent the moistur e-stress areas for which image 



classification techniques were also presented in this 
chapter. 

Application of the above techniques to limestone 
region in the present work is detailed in the succeeding 

pages . 

2.4.3 Image Classification 

2 . 4 . 3 . 1 General 

The terrain has been classified based on the image 
analysis. In the present study, three channel data has been 
used for interactive classification to make a thematic map of 
the terrain. All the four channel data have been used to 
classify the moisture-stress areas and produce a line 
printer map based on the classification. In the present work, 
supervised classification has been used. Spectral signatures 
have been studied in detail for the limestones and 
morphological features of limestones. 

2. 4. 3. 2 Supervised Classification 

It involves user in selecting some pixels from each 
possible class. From these, the classification algorithm 
determines what each class "looks like”, and then assigns 
each pixel of the image to one of the classes. In the present 
study Bayesian maximum likelihood classifier has been used. 
Bayesian Maximum Likelihood Classifier 

It is a well developed method from statistical 
decision theory which performs a decision-making function on 
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a data vector by assifinin£ it to one of a given set of 
possible classes fSwain, 1972) . The data vector can be derived 
from any set of mesurements and in the case of mult ispectral 
scanner data, the measurements are the levels of reflected or 
emitted electromagnetic energy. 

Estimation of probability with this technique 
involves two weighted factors. The a priori probability is to 
be determined first followed by a weight associated with the 
cost of misclassi f ication being applied for each class. 

The mathematical formulation involves several steps. 
The probability function is assumed and the parameters of 

that distribution are estimated. Let be random 

variables where x^ is the noise measurement of the feature. 
Let PCx/Gj^) be the conditional probability density function 
of class i and FfG^) is a priori probability of class G^. The 
task of the classifier is to assign the input sample such 
that the probability of misrecognition is minimized. 

The Bayes’ decision is that 
X s Gj^ 

that is if there are G groups then Bayes’ rule is to assign 
the feature to group i where : 

PCGi/x) > PC^j/x) for all j i ....(2.36) 

Although it is easy ’in principle' to obtain PCx/G^), the 
probability that we are interested in is P(Gi/x). Fortunately 
connection between the two, known as 
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there is a 


’ Bayes ’ 



theorm’ : 


PCx/Gi) P(Gi) 

P(Gi/x) = (2.37) 

E PCx/Gi) PCGi) 

All the items on the right hand side of the equation 
can be found by sampling. The conditional probability 
PCx/Gj^), the probability of a particular set of measurements 
’x’ given that the feature comes from group i, P(Gj[) is the a 
priori probability. 

Putting the Bayes' theorm into Bayes’ rule gives : 
assign to group i if 

PCx/Gi) P(Gi) PCx/Gj) P(Gj) 

> (2.38) 

E P(x/Gk) PCGk) E PCx/Gk) PCGk) 
all k all k 

for all j ^i 

or, as the bottom lines on both sides of the inequality are 
the same, they can be cancelled out to give i, if 

P(x/Gi) PCGi) > PCx/Gj) PCGj) for all J^i 

...(2.39) 

This is the form of Bayes’ rule. The following steps are 
neccessary to classify an image: 

(1) Training areas for each class are selected. These areas 
are known (or assumed) to contain numberof pixels in the 
class. In general the more training data per class, the 
better the accuracy, since it should give a truer 
representation of the class. 

(2) The most common choice of functional representation of 
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a Gaussian distribution in which for 


one 


P(x/Gi) is 

dimensional case, the parameters to be estimated are mean and 
standard deviation of the distribution. The general form of a 
one- dimensional Gaussian with mean m, standard deviation a, 
and normalized to have area is : 

^ C2.40) 

In the present case, it can be expressed in the form 


PCx/Gi) = e (2.41) 

'FYira j[ 

Uhere mj; and cr^ are the mean and standard "cieviation of the 

i-th class. For a multivariate n band problem, the m^ becomes 

an 'n’ dimensional vector of means, becomes an nXn 

n / 2 

covariance matrix and the scaling factor becomes (Zir) 

To perform the classification of an image using the 
parametric forms, the trainig data are used to compute ’n’ 
dimensional vector mean and nXn covariance matrix E^ for each 
class. However if within each group the variables that make 
up the mesurement vector x have a multivariate normal 
distribution, then the form of PCx/Gj^) is known. Hence 

PCx/Gj) = yT- [-HCx-mi)^ E|^ (x-mi)] .(2.42) 

C2tt)^^^|E P 
i 

In this case, estimating PCx/G^) comes down to estimating two 
paramet ers lor each group ^i ? the group mean vector and Ej_ 
the group covariance matrix- 
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Using the normal form of PCx/Gi) in Bayes’ rule gives: 


assign x to Gj^ if 


PCGi) 


exp[ -1 /2 Cx-mj_ Ej 


(x-mi)] > 


T -1 

TTn i,e5cp[-l/2(x-mj) E* (x-m;)] 

C2ti) ^^|E p -> J 

j 


for all j 


. . . . (2.43) 

Taking the natural logCln) on both sides 

of the 

inequal i ty 

gives: 



assign x to Gj[ if 



-1/2 ln(2Ti) -1/2 InjE | -l/2Cx-m )*^ E[^ 

(x-m ) + 

ln(P(G )) > 

i i 

i 

i 

-1/2 lnC27i) -1/2 InjE | -l/2Cx-m )^ E^^ 

J j ^ 

(x-m ) + 

ln(P(G )) 

j 

j 

for all j ^i 


(2.44) 


Cancelling both the terms that are common to both sides gives 

--/zlnjE I -l/2Cx-m Cx-m ) + InCPfG )) > 

i i i i 

-)^ln|E I -l/2Cx-m E^^ (x-m ) + ln(P(G )) 

j j j j 

for all j jfi ....(2.45) 


This expression has one difficulty in that the cancellation 
of all the common terms has resulted in a quantity that is 
negative. This can be corrected by multiplying both sides of 
the inequality by -1 and changing the greater than sign to 
less than sign giving: 
assign x to Gj[ if 4 



InCPCG )) < 
i 


InjE I )' Ej (x-m ) - 

i i . i 

-ln|E I Ej ^ (x-m ) - lnCP(G )) 

j j J j 

for all j fi C2.46D 

Although matrix operations are involved in working out the 
above expressions, the final result is number for which each 
group depends on x. That is 

PCG /X) = ln|E |^^2Cx-m (x-m )-lnCPCG )) ( 2.47) 

i i i ^ i i 

This is the form we use normally in the computational work, 
and in this case, p is assigned to the class for which 

P(Gj^/x) is minimum. 

Numerical Separability of the Classes 

A good classification can be obtained only if the 
classes are numerically separated. If there is extensive 
overlap in their distributions, even a high likelihood value 
for a pixel does not give us confidence in its correct 

classification. So we need some measure of class separability 
in order to assess the accuracy we can expect from a 

classification. 

Divergence is one of the first measures of 
statistical separability used in the machine processing of 
remote sensor data, and it is still widely used as a method 
of feature selection. It addresses the basic problem of 
diciding what is the best q-band subset for use in the 

t t$ 



supervised classification process. The number of 
combinations, c, of n bands taken q at a time is 


n n I 

C = 

q q! (n - q) ! 

The degree of divergence or 

classes c and d, divergence^ci . 

formula - 


(2.48) 

"separability” between two 
is computed according to the 


Diverged = 1/2 (V^ - Vd) (Vd"^ - V^"^) + 

0.5 (Vc"^ + Vd“^) (He - Md) CKc " ^d)^ 

(2.49) 

Where Tj<[.] is the trace of a matrix (i.e, the sum of 
diagonal elements, and Vd are the covariance matrix for 

two classes, c and d, under investigation, and Hd are mean 

vectors for classes c and d. 

If there are more than two classes, the average 
divergence, Di verga^.y.g , is computed over all possible pairs of 
classes c and d while holding the subset of bands, q 
constant, as follows : 



1 

m- 1 

m 


= 

II 

Z 

E Diverged 

(2.50) 


c 

c = l 

d = c+l 



It is possible to take the divergence logic one step further 
and the computed transformed divergence, expressed as 

(1 - Diverged) 

Diverge^d “ 2000 [ ( 1 - exp — ) 3 ..(2.51) 
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This statistics gives an exponentially decreasing weight to 
increasing distance between the class. It also scales the 
divergence value to lie between 0 and 2000. 

A transformed divergence value of 2000 suggests an 
excellent class separation. Above 1900 provides good 
separation, while below 1700 is a poor separation For 
divergence values less than 1500, we get spectrally similar 
classes. "5^ 
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CHAPTER 3 


DIGITAL IMAGE PROCESSING FOR LITHOLOGICAL IDENTIFICATION 
3 . 1 GENERAL 

The principal objective of enhancement techniques is 
to process a given image so that the result is more suitable 
than the original image for specific application. In the 
present work, spatial domain approach based on gray level 
mapping has been adopted. 

Enhancement procedures are carried out through 
pseudocolour concepts with superimposition procedure as well 
as through contrast stretching, histogram modification and 
edge enhancement techniques. Uhile former procedures have 
enabled distinct demarcation of lithological variations, the 
latter have enhanced the structural features such as fracture 
zones, lineaments as well as the associated sinkhole 
features. The interpretations have been correlated with 
hydrogeological field data to confirm the compatibility of 
the techniques in detecting these features. Software has been 
developed interactively for all these techniques by the 
author. Few areas were selected, for enhancement, which 
apparently represented dense vegetation under moisture-stress 
conditions, thought to have been induced by sinkhole 
collapsed formation. Using enhancement techniques, various 



litholofiical units have been delineated for further 
hydrogeologic evaluation. Extensive field observations were 
undertaken as ground checks. 

3.2 LITHOLOGIC IDENTIFICATION 

For this purpose, enhancement procedures such as 
contrast stretching, pseudocolouring and superimposition have 
been used. Compactness of the limestone, its tonal contrast, 
erosional features, and the drainage pattern are the 
parameters which aid in the terrain differentiation. Based on 
these parameters, it is possible to demarcate the lithologic 
units. Distinct lithological differences within the 
limestones have been recognized in the present study. In 
addition to these variations in the limestones, structural 
features of hydrogeological significance have also been 
delineated through enhancement techniques. The study area can 
be broadly divided into three lithological units : 

(i) Limestones. 

(ii) Shale. 

(iii) Alluvial tract. 

To distinguish them distinctly the enhancement 
techniques used are ratioing, histogramming , histogram equa- 
lization, contrast stretching and edge enhancement. 

3.2.1 Limestone 

This unit, designated as the Chandi formation within 
Raipur group, occurs in the central part of the Chhattisgarh 
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basin. The Chandi formation has been subdivided into Nipania, 
Pendri and Newari limestone members. Distinct differences 
within these members have been recognized on the basis of the 
above mentioned enhancement procedures. For this purpose, 
pseudocoloured superimposed method and ratioing and contrast 
enhancement techniques have been used. The Nipania member is 
characterized by its light tonal contrast (indicated by ’N’in 
Figs. 3. 8 and 3.11) in pseudocoloured photographs. It could be 
noticed that the oval features such as sinkholes are very 
scant in the limestone due to its compactness. This member is 
located on the N-U side of the Jamunia river between 
Hathbandh and Bhatapara areas. On the other hand, on the 
eastern side of the same river in the region south of Tilda, 
the Pendri and Newari members were identified by their dark 
tonal contrast and soluble nature of the limestone. For a 
distinct identification, ratioing and contrast stretching 
methods have been utilized. 

In the ratioing process, two spatially registered 
spectral images are divided pixel by pixel (Siegal, 1980). 
Ratioing can be thought of enhancing minor differences 
between formations (between Pendri member and Nipania member 
in the present case). The Pendri member can be clearly 
identified by its dark gray tonal contrast with the dominant 
occurrence of oval and irregular features (indicating ’P’in 
Figs. 3.9 and 3.12) revealing the soluble nature of the 
member. 
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To confirm that all the dark irregular features in 
not represent the sinkholes, contrast enhancement 
has been done on the ratioed (2/4) picture. For obtaining the 
contrast enhancement, the probability density function of the 
reflectance values of the area of interest should be known. 
For this purpose, a histogram was plotted for the area 
(indicated by ’I’ in Fig. 3. 14). It has been observed that 
most of the reflectance values in the histogram were 
distributed between 1.0 to 2.5. These values were later 
subjected to stretching for the entire medium to improve the 
contrast of the area. Through this enhancement, the actual 
sinkholes were identified clearly by their dark tone 

(Fig.S.l?*). The gray tone linear feature (indicated by ’LV’ 
in Figs. 3. 14 and 3.1'?') is a recent active zone along which 
the solution action is prominent. The rest of the areas with 
light gray tone in the vicinity reveal dense vegetation 
along the fracture zone indicating the proximity of 
groundwater. The exposed surface of the Pendri member has 
been identified by its light tonal contrast (shown by ’P’ in 
Figs. 3. 11 and 3.12). Compared to the Newari member, this 
member has greater solubility. The Jamunia river follows the 
contact of Nipania and Pendri members. This area was later 
subjected to selective stretching to enhance the subdued 
caverns under thick vegetational growth (indicated by ’CV’in 
Fig.3.2.0.j. The red patches indicate the vegetational growth 
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and distribution pointinfi to the proximity of the 
£roundwat er . 

The third member has been identified by its light 
and gray coloured irregular shape (indicated by ’NE’in 
Figs. 3. 9, 3.12, and 3.16). The river Banjari has its initial 
course in the Newari member and the remaining part is along 
the contact of Pendri and Newari members. The Newari member 
(indicated by ’NE’ in Fig. 3. 9) is also a soluble limestone. 
In this member, several collapsed zones and sinkholes have 
been identified. To distinguish this member clearly, ratioing 
and contrast stretching have been used. The exposed portion 
of this member is restricted, which can be identified by its 
light tonal contrast with irregular shape (Figs. 3. 9 and 
3.12). The ratio of bands 2 and 4 was used to detect the 
sinkhole and collapsed zone areas. They have been located 
with irregular and oval shape features in the northern part 
of the area. A histogram has been plotted for these areas to 
know the probability distribution function of the reflectance 
values and it has been observed that most of the values have 
been in the 2 to 3.5 range. This range has later been 
stretched to entire medium to enhance clearly the collapsed 
zones and sinkholes. These were represented by ’CZ’ in 
Figs. 3. 15 and 3.18. Edge enhancement has enabled 
identification of the actual collapsed zones filled with 
sediments and with dense vegetation zone from within the dark 
tonal contrast zone, the entire portion which apparently 
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looks other wise like a collapsed zone. The oval collapsed 
zone with a storage tank is a recent one located on eastern 
border of the Newari member and has been identified by its 
gray tonal contrast. Uithin this limestone member, it has 
been observed in the present study that the groundwater is 
under unconfined conditions in the sediment filled area of 
the collapsed zone, while in the deeper regions, where 
intercalations of Gunderdehi shale occurs, the confined 
conditions exist. The surface of the limestone terrain is 
characterized by discontinuous drainage system in the 
vicinity of the oval features. 

3.2.2 Shal es 

Two shale members, the Raipur and Tarenga, have been 
identified in the imageries and were later confirmed by field 
checks. The Raipur shales are recognized in the contrast 
stretched window (indicated by ’Sh’ in Figs . 3 . 9 , 3 . 12 and 
3-18) along the Jamunia river. These shales, interbedded with 
Pendri limestones and Nipania limestones, are poorly visible 
in bands 2 and 4 due to paucity of their exposure and 
extensive vegetation. Checks around Hurpar and Diggi villages 
have enabled the confirmation of these shales (Figs. 3. 2*^ and 
3. So). To enhance these areas to the maximum extent, 
selective stretching procedure was applied by choosing the 
suitable stretch parameters. The Raipur shales are 
characterized by the dark tonal contrast (indicated by ’ Sh 



in Fift.3.18) due to their f ine~grained nature with presence 
of moisture. 

The Tarenga shale occurs in the region north-vest 
and north of Ti Ida-Bhatapara line (indicated by 'Tsh' in 
Figs. 3. 2 and 3.8) . Although, in general this shale is also 

of fine texture, it exhibits relatively lesser tonal contrast 
due to the absence of moisture as compared to the Raipur 
shales. 

3.2.3 Alluvium 

The major alluvial tract of the terrain has been 
identified in the north-western part of the image near Rajpur 
where the river Jamunia has its confluence with Seonath river 
(indicated by ’Al’ in Figs.3.4, 3.7, 3.10 and 3.13). To 
demarcate the area clearly, it has been subjected to ratioing 
and selective stretching. The fine-textured gray coloured 
tonal contrast, as indicated by ’Al’ in Figs. 3. 16 and 3.19 
aids in the interpretation of this tract. For the alluvial 
tract, reflectance is a sensitive indicator of the terrain 
along with the vegetational growth. Variation in reflectance 
values is related to the compositional variations of the sub- 
surface in the alluvial regions. The alluvial patches are 
potential sources for groundwater in the region and are 
extensively developed. 

Edge enhancement procedure has been utilized to 
demarcate the structural features such as fracture zones and 
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lineaments and to demarcate perfect boundaries of sinkholes 
and collapsed zones. It is mainly based on high~pass 
filtering through which sudden reflectance changes can be 
delineated and edged. This would facilitate the demarcation 
of various structural features. The software for the edge 
enhancement has been developed by the author during the 
present study. 

Laplacian 3X3 and 5X5 filters have been used in 
the present work for a convolutional operation. The edge of a 
sinkhole or a collapsed zone represents some kind of 
irregularity with respect to the surrounding pixels. Due to 
this change, the sharp edges of the limestones (Fig. 3. 21) 
were delineated. The measure of the magnitude of the 
irregularity at a specific point, is the difference between 
the actual value for that point minus an interpolated value 
based on the remaining pixels. The actual sharp edges of the 
sinkholes and caverns were edged and assigned a black colour 
to differentiate from the surrounding areas (indicated by ’CZ 
and SK ’ in Fig. 3.21) within the region covering Tilda, 
Hathbandh and Nawapara villages. This area has been inferred 
as active solution action zone which was later confirmed by 
ground checks . 

The present investigations have revealed the maximum 
occurrence of the collapse features .in the area between 
Pendri and Arjuni villages. All the collapsed and filled 
zones and sinkholes were edged due to sharp differences in 



reflectance values and were assigned a black colour for these 
areas (Fig. 3. 21) to represent the structurally —controlled 
moisture zones. Dense vegetational growth has been identified 
along the fracture zones indicating the proximity of the 
groundwater (demarcated by ’V'). Vegetation has been assigned 
red colour for easy identification. The shaded red areas 
indicate the slopes of the limestone areas with vegetation. 
The major fracture zone with NE-SU trend along which the 
Jamunia river flows, has been edged due to sharp change in 
reflectance and was identified by its linear nature 
(indicated by L-L in Figs. 3. 23 . ) . To obtain clarity 

of information about the boundaries of limestone edges, the 
edge enhancement procedure is performed on original band-4 
windowed data. The white coloured curved edges represent the 
boundari es (B) of the limestones and straight lines represent 
lineaments and fracture zones (indicated by L-L in Fig. 3. 2$)- 
These weak zones have been confirmed through field checks. 
Pumptests in the present work have revealed that such areas 
are potential areas of groundwater yield with high 
transmissivity values. 
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ig.3.1 Enhanced window indicating the fractured limestone 



ig.3.2 Enhanced window of band 2 (Bhatapara-Tilda region) 
C = Upland tract; Tsh= Tarenga Shale 

— 8^5 ^ 



Fig. 3. 3 Enhanced window of band 2 (Ar juni-Pendri region) 


Upland tract; M = Moisture zone 


Enhanced window of band 2 (Rajpur region) 


Fig. 3. 4 
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Fig. 3. 5 Enhanced window of band 4 (Bhatapara-Tilda region) 

C = Upland tract; H = Hoisture zone; 

1 

R1 = Railway line 



Lg. 3.6 Enhanced window of band 4 (Ar juni-Pendri region) 

C = Upland tract; V = Vegetation; Sk = Sinkholes 
2 

M = 


Moisture zone 
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Pseudocoloured and superimposed section of band 2 

CAr juni-Pendri region). C = Upland tract; 

2 

M = Moisture zone; WE = Newari limestone; 

P = Pendri limestone; Cz = Collapsed zone; U = Water 



Pseudocoloured and superimposed section of band 2 

(Rajpur region) . C = Upland tract; A1 = Alluvial tract 

1 





11 Pseudocoloured and superimposed section of band 4 


(Bhatapara - Tilda region). M = Moisture zone; 

N = Nipania limestone; Sk = Sinkhole; R1 = Railway line 
P = Pendri limestone 



g.3.12 Pseudocoloured and superimposed section of band 4 
(Ar juni-Pendri region) . M = Moisture zone ; 

NE = Newari limestone; P = Pendri limestone; 

= Raipur Shale;yjpz = Gol lapsed zone; 



3.13 Pseudocoloured and superimposed section of band 4 
CRajpur region). A1 = Alluvial tract; 

N = Nipania limestone 



i.l4 Band ratio (2/4) with histogram (Bhatapara-Tilda 

region). 11= Moisture zone; N = Nipania limestone; 
Lv = Linear feature with vegetation 






Fig. 3. 15 Band ratio (2/4) with histogram (Ar juni-Pendri region) 


WE = Newari limestone; P = Pendri limestone 


UT = Surface water tank within the collapsed zonefCz) 


Fig. 3. 16 Band ratio (2/4) with histogram (Rajpur region) 


A1 = Alluvial tract 



Fig. 3. 17 Linear stretching with histogram (Bhatapara-Tilda 


region). M = Nipania limestone; Lv = linear feature 


NE = Newari limestone; P = Pendri limestone 


Sk = sinkhole 





Fig .3.19 


Linearly stretched and enhanced part of the Jantunia 
river (Rajpur region). Al= Alluvial tract 



Fig. 3. 20 Selective stretching of part of the Jamunia river 
CBhatapara-Tilda region); Cy = cavernous zone 




Fie. 3. 21 Edge enhancement section (Bhatapara-Tilda region) 
Cz = Collapsed zone; Sk = Sinkhole; 


M = Moisture-stress zone 


Fig. 3. 22 





Edge enhancement section (Arjuni - Pendri region). 

Red patches indicates dense vegetation (V) in fracture 
zones and dark patches indicate moisture zones (M) and 
sinkholes (Sk) and Collapsed zonefCz) 


Fig. 3. 23 


Edge enhancement section (Arjuni - Pend; 
B = Limestone edged boundary; L-L = Lint 



Fig. 3. 24 Stretched and enhanced part of the Jamunia river 
(Ra jpur-Ar juni region). F F = Fracture zone; 

Cz = Collapaed zone 



Fig. 3. 2 5 S^tretched and enhanced part of the Jaiaunia river 
CBhatapara-Ti Ida region). M = Moisture zone; 

F--F = Fracture zone 
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Fie.3.26 Thematic map representing various land cover types 


in the Study area 


Fig. 
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Histogram equalization of the study area 



26 Hi^'tOi^ram equal iisat ion with vegetational enhancement. 
All red patch areas represent vegetation 







ur Shale iia the course of Jamunia river near Murpar 







CHAPTER 4 


digital image processing for morphology. 

DRAINAGE AND LAND USE 

4.1 GENERAL 

To study the morphology and drainage features, 
windowed noise-free images and pseudocoloured superimposed 
procedures were used. The terrain in the central portion of 
.{•he Chhattisgarh basin has distinct and characteristic 
variations in terms of sinkholes, fracture zones and 

drainage. A study on regional scale of the fractured and 
collapsed zones in the limestone formations would facilitate 
the recognition of changes in the morphological units within 
the terrain. As the area (Fig. 3.1) undertaken in present 
study being extensive and covered with vegetation over the 
limestone terrain, digital data of Landsat-4 has been chosen 
for analysis and interpretation. Different morphological 
units have been characterized by the changes in their 
reflectance values. In addition to the linear structural 
features. drainage characteristics are recognized on a 
regional scale. The processed picture data aid in the 
delineation of sinkholes and fracture zones within the 

limestone and to understand the nature of groundwater in 
terms of its proximity to the surface. In the present 
chapter, the analysis of the processed data has been carried 
out on the basis of which, the morphological and structural 



features have been demarcated. Various sinkhole and collapsed 
zones have been mapped to indicate their association with 
areas of groundwater availability. An attempt has been made 
to determine the water table depth using reflectance values 
within the collapse zones filled with sediments as also the 
moisture zones. 

Different techniques used in the present study have 
been described in the previous chapter. After checking all 
possible combinations, bands 1, 2 and 4 have been considered 
best for the present study. 

The limestone and shale contacts have been clearly 
visible in the enhanced windowed data. Based on the tonal 
contrast, various morphological features have been identified 
along with distinct boundaries. Surface water bodies have 
been mapped using band-4. Areas with dense vegetative growth 
indicative of proximity to groundwater have been mapped using 
superimposition procedure. Ground checks have been carried 
out at various places for checking the interpretation. The 
results are presented in this chapter. 

4.2 MORPHOLOGY 

Hydrogeological deciphering is possible by studying 
the morphological features with characteristic relief. Since 
f'^^.ief gives an idea about the distribution of surface water 
infiltration, recharge, discharge and also the flow of 
groundwater, study of the typical relief forms in a terrain 



provide an understanding of the groundwater conditions of 
that terrain. Morphological units have been Snapped on the 
basis of tone and erodability characteristics of the 
limestones. 

The study area has been divided on the basis of 
morphological features into 

(i) Upland tract 

(ii) Valley tract 

(iii) Moisture-stress areas 

The upland area is characterized by its flat 
topography with thin soil cover. The valley tract exhibits 
sharp cutting nature throughout the central region. The 
moisture-stress zone is situated in the eastern and north- 
eastern part of the Jamunia river. To delineate the various 
morphological features of the study area, bands 2 and 4 as 
also the enhanced noise-free photographs along with 
pseudocoloured products have been used. 

4.2.1 Upland Tract 

This tract in general has flat topography and 
exhibits variation in tonal contrast in the imageries due to 
variations in soil texture and groundwater situation. Uhen 
compared to the adjacent low lying areas, this tract is 
lighter in its tone due to its high reflectance. Uithin this 
tract, distinct differences in reflectance values enable a 
subdivision into two clear regions (designated as and C 2 
in Figs. 3. 2 to 3 . 6 and 3.8 to '3.10). The tract between 
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Batapara and Hathbandh (Cj^) could be recognized by its tone 
due to high reflectance values (ranging between 70—77) 
compared to the tract ’C 2 ’near Arjuni (ranging between 65— 
71). The tract ’C 2 ’ is evidenced with frequent intermixing of 
upland and ravinous-1 ike features and exhibits a wide 
variation in its erosional status, as in the areas around 
Arjuni and Mopar villages. For example in the tract around 
Arjuni, erosional features such as depressions and gullies 
with gently rounded cross sections (Figs. 3. 3 and 3.9) are 
distinctly visible. The fracture zones in the upland tract 
offer channels for surface water percolation to underground 
streams. The natural tanks, oval in shape, are indicative of 
the sinkholes with surface storage. Groundwater in this 
region is tapped through dug wells in fractured and jointed 
zones. The hydrogeological details of the same are presented 
in the subsequent chapter. 

4.2.2 Valley Tract 

This tract is characterized by its intense erosion 
due to solution action and associated erosional features such 
as gullies with intermittent ridges that distinctly stand out 
in relief. The ravinous gullies might have formed due to 
solution action along joints and fractures. The erosion 
action is maximum throughout the valley area. This erosion 
due to solution action continues its operation until 
groundwater is reached. Thus the depth of gullies in the 



valley tract indicates the depth of groundwater. The deep 
valleys with a dry appearance can be seen clearly in the 
pseudocoloured figures (indicated as ’VL’ in Fig. 3. 8). The 
with the Banjari river channel is characterized by a 
wide flat terrain underlain by shale. The impervious nature 
is clearly indicated from the flow of surface water (as 
designated by ’U’ in Fig. 3. 3). Field checks confirm such an 
observation. 

4.2.3 Moisture-stress Areas 

These areas have been delineated on the generated 
images on the basis of their dark tone (Figs . 3 . 2^ and 3.25!) 
reflecting the presence of finer nature of the sediments with 
proximity to groundwater. For example, the irregular areas 
with dark tones around Nawapara and Kesili villages 
(Figs. 3. 11 and 3.14 as indicated by ’M’) as well as around 
Ar juni are the moisture-stress areas which have also been 
confirmed through ground checks. Here the groundwater occurs 
in both confined and unconfined conditions. The coarser 
pediment zones that are prominent locations for good 
groundwater conditions, could be clearly indicated along the 
slopes of the upland tract. Surface water tanks and filled 
zones have been noticed on band-4 and pseudocoloured 
imageries ( indicated by ’M’ in Figs. 3. 9 and 3.12). 
Extensive vegetation in these areas reflects proximity of 
the groundwater. A number of collapsed zones have been mapped 
in the eastern part of Banjari river. One of such collapsed 
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zones with distinct oval outliers with water storage tank is 
identified and indicated in the Baloda Bazar area (indicated 
by CZ * in Figs. 3. 9 and 3.12). On the Western part of this 
ringed area, an irregular collapsed zone of a relatively 
older age has been identified with heavy vegetation 
(indicated by 'M’ in Fig. 3. 12). Another ring-shaped unit has 
been located in the northern part of the imagery which might 
have been recently formed by intermixing of sinkholes with 
vegetation (indicated by in Fig. 3. 12). Dense vegetational 
areas were enhanced through histogram equalization and false 
colouring (Fig.3.2B^) representing the areas of groundwater 
availability in the proximity. Tube wells in few of these 
localities have been operating with good yield. Field checks 
have confirmed dense vegetation in these areas. The 
hydrogeological aspects of these areas are detailed in the 
subsequent chapters . 

4.3 DRAINAGE AND STRUCTURAL CONTROL 

Drainage pattern plays an important role in 
influencing the hydrogeological situation of the limestone 
terrain. In general, drainage conforms to the regional 
structural trends. In the present work the structural control 
of the drainage pattern has been established on the basis of 
observations of the drainage network as related to the 
fractures and lineaments on the satellite imageries. The 
drainage pattern largely is of rectangular type. The Jamunia 
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river itself is flowing in a fracture zone trending NE-SU. 
Initially the river has a N-S trend which suddenly changes to 
east (Fig.S.S) following a fracture zone (indicated by ’F’ in 
the Fig. 4.1). In this fracture region, considerable amount 
water infilters, contributing to solution action. 
Consequently, caverns of various sizes have formed (indicated 
by ’S’ in Fig. 4.1) through which the water flows internally 
at the place indicated by ’Aj ’ in Fig. 4.1, and emerges again 
at ’^ 2 ’ surface. Similar feature also been evidenced 

again at another place (indicated by in the Fig. 4.1) 

infiltering through a fracture zone and subsequently emerging 
at ^ 2 - In fact, all the streams contributing to the Jamunia 
river are controlled by regular fracture patterns. The 
drainage network in general is very coarse in texture with 
considerable amount of water infiltering through fractures 
and solution channels. Near Kesili and Nawapara, caverns of 
large size and sinkholes have been identified along the river 
(indicated by an arrow in Fig.3.er) reflecting the soluble 
nature of the limestone. In the eastern and southern parts of 
the river (Figs. 3. 6 and 3.12) surface water bodies, such as 
tanks are noticed as dark isolated spots in the terrain. Some 
of the tanks were connected to the natural drainage. The 
major part of the Jamunia river course appears to be 
controlled by these .structural features. Present study 

revealed that the groundwater in these areas is under 
confined conditions. 
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4.4 IMAGE CLASSIFICATION 


Supervised classification method, Bayes maximum 
likelihood classifier, has been used to classify various land 
use in the study area. The details of the classifier have 
been elaborately dealt in chapter 2. The statistical analysis 
of the training set of the classes given in Appendix III. The 
parameters calculated for all classes in four bands are, 
mean, variance-covariance matrix, inverse of variance 
covariance matrix and correlation coefficients of all classes 
in four bands. 

Once the above parameters are calculated for a 
training set, we may take reflectance values in all four 
bands of any unknown feature based on the logic described in 
chapter 2. 

Initially training set itself is subjected to the 
classifier giving 40 samples for each class. Classes 1,3, and 
5 have been classified 100 percent except class 2 which is 
97.5 percent. The confusion matrix has given in the results 
showing the accuracy of the above mentioned classification. 
Infact before selecting the classes , the data has been 
checked with the numerical separation technique, which gives 
the idea about the overlapping. A divergence/transformed 
divergence matrix has been calculated to represent the 
separability within the classes. As can be seen from the 
results (Appendix III) all the classes mentioned above were 
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easily separable (e.fi. classes 1,3,4 and 5 are almost near to 
2000, but the class 2 has around 1850 which may be due to 
one misclassi f icat ion shown in confusion matrix). Later the 
entire study area has been classified based on the same lo£ic 
and a line printer map (Fig.4.2) has been produced to 
represent the classes. Spectral reflectance curves (Figs. 4. 3 
and 4.4) have also been plotted for various land use of the 


study area. 



CHAPTER 5 

GROUNDWATER SITUATION AND AQUIFER CHARACTERISTICS 


5.1 GENERAL 

The carbonates in the central part of Chhattisgarh 
basin belong to fourth cycle (Schnitzer , 1 96 9 ) of the 
carbonate-shale sequence. These limestone formations have 
been demarcated as three distinct members namely Nipania, 
Pendri, and Newari members . A comprehensive hydrogeologic 
investigation of these limestone aquifers includes not only 
the distribution, structural features favourable for the 
recharge and storage of groundwater but also their 
performance characteristics as aquifers. The potential 
aquifer zones within the members together with all the 
associated structural features have been delineated on the 
basis of the remote sensing analysis as detailed in the 
previous chapters. Investigations have been carried out on 
groundwater situation and the seasonal fluctuations of water 
levels in the limestone aquifers- Systematic pumptests have 
been conducted in selected areas, delineated through 
enhancement procedures and the aquifer characteristics have 
been worked out. A brief study on the chemical quality of the 
water has also been undertaken to provide complete picture of 
the groundwater hydrogeology of the study area. The present 
chapter contains the details on all these aspects. 



5.2 GROUNDUATER SITUATION 


Primary porosity plays an important role in makin£ 
the limestone formations as good aquifers. If the formation 
has little primary permeability with sluggish water 
circulation, it will be a poor aquifer. On the other hand if 
the rocks are in the zone of active circulation of water, 
they are usually cavernous and are likely to be good 
aquifers. 

The carbonate rocks in the present area are of 
Raipur group, comprising of Chandi formation. The Chandi 
formation has been classified into three limestone members- 
(1) Nipania (2) Pendri and (3) Newari. Of the three members, 
the Newari member is less soluble being massive and devoid of 
joints. The recharge area for this member is in the eastern 
side where it has a gradational contact with Gunderdehi 
shale. Uithin this member in the area around Baloda Bazar on 
the northern side, collapsed zones with storage tanks were 
noticed, and these are formed by underground channeling from 
gradational contact. A collapsed zone with an irregular shape 
has been identified near Arjuni. This collapsed zone was 
later filled by recent alluvium which forms a potential 
source for groundwater occurrence within the member. The 
Pendri member is a soluble limestone. This formation is 
extensively jointed, fractured and has stromatolites. Shale 
intercalations are noticed towards the top surface of this 
member. Extensive solubility in this limestone on a wide 


U4 



scale has been noticed along the fracture Eones as well as 
along the contact of the stromatolites connected by carbonate 
mud. Sinkholes, potholes and collapsed zones are more 
prominent in this zone. Sinkhole distribution appears to have 
been controlled by weak zones associated with fractures. 
Surface drainage contributes to a substantial part of the 
recharge through fractures in the western and north-western 
region. As this member is more cavernous and has fractures 
and lineaments, it acts as a storage reservoir with high 
yield. 

The Nipania limestone is coarse to medium-grained 
with mottled appearance. Fractures and lineaments are 
scantily distributed. This member, at a relatively higher 
elevation compared to the other members, has its recharge 
from monsoon. On the eastern side of the member, fracture 
zone helps in the recharge which contributes in the 
groundwater storage across the boundary to the Pendri 
limestone member. Sinkholes and collapsed zones are 
conspicuously absent in this member. 

The limestone members in general are characterized 
by secondary porosity and permeability controlled by 
fractures and solution channels. Since dissolution is 
localized along the weak zones, such porosity and 
permeability is highly directional and strongly anisotropic. 
Using remote sensing techniques, the fracture traces and 
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lineaments were delineated which act as a guide for locating 
the directional zones of porosity and permeability and hence 
could be used as a guide for locating the high capacity wells 
in these limestone members. 

The main source of recharge to the groundwater in 
the study area is through precipitation and surface run-off. 
The surface run-off that approaches the sinkholes and other 
openings and channels in bed rock quickly replenishes the 
groundwater. It has been observed that the surface drainage 
as it enters the cavernous areas infilters to the subsurface 
as evidenced from sudden depletion of water in the channel 
(Fig. 4. 1-in the span from to A 2 ). The water again emerges 
into the channel as soon as the cavernous region is passed. 
This is clearly evidenced in Fig. 3. 9 where the channel on the 
northwestern corner has water (indicated by the lighter tone) 
while in the cavernous region (Fig. 3.8) same channel in its 
upstream portion has no water, (as evidenced from its dark 
tone in the central part of the picture). The lithological 
and structural variations also control the depth to water 
level in the area. The movement of groundwater in carbonate 
rocks is complicated by variation of permeability within the 
limestones. To demarcate the recharge and discharge areas in 
the study area, water level maps (Figs. 5.1 and 5.2) have been 
prepared for both pre— monsoon and post— monsoon periods. The 
arrows in these maps indicate the possible flow direction of 
the groundwater. The flow of water has gradient from west and 
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northwest and northern regions sloping toward the central 
region where the Pendri member is located. Water level 
fluctuation map (Fig. 5.3) prepared on the basis of water 
level data in the pre-monsoon and post-monsoon periods has 
revealed maximum fluctuations wherever the sinkholes, caverns 
and collapsed zones are present. For example the fluctuations 
are maximum around Bitkuli, Arjuni and Jhipania and also 
within a substantial part of the central region where soluble 
and cavernous Pendri member is present. In general, the 
groundwater topography indicates a flow in the direction of 
river flow. However, there are extended patches (Figs. 5.1 
and 5.2) where the groundwater slope appears to be very small 
indicating the likelihood of storage in caverns. Due to the 
presence of caverns acting as sinks, we cannot expect one to 
one correspondence between the groundwater topography and the 
river flow direction. 

Concentrations in the groundwater samples from 
different locations together with the values of conductivity 
and pH are given in Appendix-IV. The general increase in the 
hardness in the direction of groundwater flow noticed can be 
attributed to the increased duration of the contact of water 
in caverns and filled zones within limestone member. In the 
recharge areas, the higher concentration of CO 2 and low 
dissolved solids content promotes dissolution of calcite and 
dolomite. The water in the study area in general can be 
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classified as weakly alkaline and soft to moderately hard. 
Existence of numerous ponds in the area appears to have a 
pronounced affect on chemical quality of groundwater in the 
wells. It has been observed that locations close to ponds 
show high HCO 3 and corresponding hardness as evidenced in 
Arjuni, Bhatapara and Bothidih villages. 

5.3 METHODOLOGY OF PUMP TESTS 

5.3.1 General 

In the present work, aquifer evaluation has been 
carried out through pumptests on tube wells penetrating 
limestone aquifers in selected areas delineated by the author 
through remote sensing analysis, and time-drawdown as well as 
time-recovery data have been collected. As there is no 
observation well facility, all the pump tests were carried 
out in the abstraction well itself and the data has been 
analyzed by Theis, Jacob and Ualton methods. 

5.3.2 Analysis and Interpretation of Aquifer Tests 

The hydraulic functions of an aquifer are to store 
and transmit water. The ability to store water is expressed 
by the storativity, S, defined as the volume of water that an 
aquifer releases from or takes into storage per unit surface 
area per unit change in the component of head normal to that 
surface. For water table aquifers, S is equivalent to 
specific yield of material, and its values range between 0.01 
and 0 . 35 ; values for artesian aquifers range from 10 to 
10 "^ CHalton, 1964). 
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The ability to transmit water is 


expressed by 

transmissivity, T, defined as the rate of flow of water, at 
the prevailing temperature, through an unit width of the 
aquifer extending over the entire saturated thickness of the 
aquifer, under hydraulic gradient of 100 percent, which 
means 1 m drop in head over 1 m of flow distance. 
Transmissivi ty is numerically equal to the product of 
permeability of the material and the saturated thickness of 
the aquifer. In addition to the transmissi vi ty , the 
storativity has also been calculated, using Ualton's specific 
capacity equation, for the wells in the carbonate terrain. 

All the test data analyzed have been obtained on 
the tube wells penetrating limestone and shale aquifers. 
Discharge from the wells has been measured with 90® V-notch 
or with piezometer depending upon the convenience and the 
facilities available in the field at individual well sites. 
The water level fluctuations during the drawdown and recovery 
periods have been monitored in the wells. 

In the tests conducted, the static water level has 
been measured in the well in each case prior to the 
commencement of pumping operation. The drawdown measurements 
are made since the start of pumping at different time 
intervals during the pumping. These intervals have been 
chosen depending upon the pattern of variation in the 

well. Pumping has ; been carried out 


drawdown in the 


for 



duration ranging from 250 min. to 6 hours, depending upon the 
availability of the wells. The recovery measurements were 
carried out on the cessation of pumping and the time 

intervals in recovery measurements are chosen to reflect the 
true pattern of variation of water levels during recovery 
■period. Time — drawdown data have been interpreted using 

Jacob’s straight line method and the recovery data by Theis 

(1935) method. To determine the storativity, Ualton’s 

specific capacity equation has been utilized. 

5.3.3 Theis Non-equilibrium Formula 

Theis (1935) derived an equation, the non- 
equilibrium or Theis formula, which describes spatio 
temporal effects of pumping a well in an aquifer. Uhen a well 
penetrating an aquifer is pumped at a constant rate, the 
influence of discharge extends outward with time. In the case 
of confined aquifer, the rate of decline of head multiplied 
by storage coefficient and summed over the area of influence 
equals the discharge. As the water is drawn from the storage, 
the head in an aquifer continues to decline for sufficient 
period of time and thus it is not possible for steady state 
flow to be attained. 

The Theis non-equilibrium equation is written as 

m — tj: 
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and consequently 
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S = 


(5.2) 


s = drawdown, in m, measured in a piezometer at a distance 
’r’ from the abstraction well, 

3 

Q = constant well discharge in m /day, 

S = storativity, dimensionless, 

T = transmissivity of the aquifer in m*/day, 
and 

t = time in days since pumping started. 

The value of the integral expression in equation 
5.1, also written as U(u) and read as ”Uell function of u,” 
is given by an infinite series. 


**> -u 

e 


du 


u u 


UCu) = - 0.577216 - log u + u - + 

e 2.2! 


(5-3) 


Values of ’U(u)’ for values of ’u’ as tabulated in Uenzel 
(1942), are given by Ferris and others (1962). I f the values 
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for UCu), Q and s are known, value of ’T’ can be calculated 
from equation 5.1. However, to know U(u), both ’s’ and ’T’ 
roust be known, and in the equation 5.1, 't* occurs both in 
the argument of the function and divisor of the exponential 
integral . Hence ’T’ and ’S’ can not be determined directly 
from this equation. However, a graphical method of 
superposition, popularly known as "type curve method”, 
devised by Theis and described by Jacob (1940) is used to 
determine ’T’ and ’S’. 

From the equation 5.1 and 5.2 it follows that 

Q 

log s = log UCu) + log (5.4) 

4 Ti T 

1 r* S 

log t = log + log (5.5) 

u 47r T 

Equation (5.5) has been arranged for having 1/u and t by 
Ualton instead of u and r*/t of Theis. For a constant 
discharge Q, the equations 5.4 and 5.5 show that U(u) is a 
function of (1/u) in the same way that ’s’ is a function of 
t. Therefore, it is possible to superimpose the data curve (s 
vs t plot) on double log sheet of the same modulus as type 
curve on the Ualton’s type curve (U(u) vs 1/u plot on double 
log sheet), keeping the corresponding axes of the graphs 
parallel to each other and in such a way that the data best 
fit the type curve. The match point coordinates U(u), 1/u, s 

and t are read from the respective graphs, which are used in 
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equations 5.1 and 5.2 to estimate the aquifer parameters T 
and S . 

5.3.4 Jacob’s Method 

Jacob (1950) modified equation 6.1 to devise a 
straight line graphical method for the estimation of 
transmissivity and storativity. This method is based on the 
fact that when u becomes small (u<0.01) the sum of terms 
beyond three terms (Eqn. 5.3) is not significant and the 
equation 5.1 reduces to - 


2.3 Q 2.25 T t 

s = log r (5.6) 

4 Tt T r S 

This shows that a plot of drawdown against the 
logarithm of time after pumping started describes a straight 
line. The slope of the line is equal to 2.3 Q/4'nT. For a 
pumped well, if the change in drawdown (s=S 2 -sj) is taken 
over one log cycle of time (logiQt 2 /t 3 ^ = l ) , then the slope 
becomes 


and 


when s 
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0, the equation 5.6 becomes 
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which gives 


2.25 T t 


o - — (5.11) 

r* 

where S,T, and r are previously defined and t is the time 
intercept; that is, the’ time at which drawdown (s) is zero. 
Equation 5.11 is used only with data obtained from 
observation well. 


when u = 0.01 , equation 5.1 gives 


r^ S 


t = 


0.01 X T 


(5.12) 


or 

10^ X r^ S 

t = (5.13) 

T 

Equation 5.12 gives the time, in minutes, that must elapse 
before straight line method can be applied to pumping test 
data. This shows that t increases as r and S increase, and 
decreases with increasing T. For a constant value of S/T, t 
increases with r and for an abstraction well, where r is the 
radius of well, t is the smallest. 

Using equation 5.12, the periods of time that must 
elapse before the straight line method can be applied to a 
pumped well of 7.6 Cm radius, for several assumed values of S 
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and T were calculated ( A . -fhis shows that except for 
five values, which are for combinations of T and S, ' rarely 
observed in the field, strai£ht line method is applicable 
after a time period rangin£ from few seconds to 16 minutes of 
pumping. This indicates that the time-drawdown data from the 
pumped wells could be analyzed, using straight line method, 
to calculate the values of transmissivity. 

5.3.5 Theis Recovery Method: 

The recovery of water level in a well since the 
pumping is stopped is expressed in terms of residual drawdown 
s, which is the difference between the static water level 
and the actual water level measured at a certain time t since 
pumping stopped. The recovery method has an advantage in 
that, due to the constant rate of recharge ’Q’ being equal to 
the mean rate of discharge ’Q’ during pumping, the drawdown 
variations resulting from slight differences in the rate of 
discharge do not occur during recovery (Kruseman and Ridder, 
1970). The residual drawdown, s’, during the recovery period 
according to Theis(1935), given by: 

t 

Q 4 T t 4 T t 

s’ = In - In r- (5.14) 

4 Ti T r* S r* S 

where 

9 ' 

s = residual drawdown in m. 


r 


distance 


in m from pumped well to observation well 



(if observations are recorded in the pumped well 
itself, r = ry = effective radius of the pumped well) 

f 

S - storativity during recovery, dimensionless, 

£ - storativity during pumping, dimensionl ess , 

2 

T = transmissivity in m /day, 

t = time in days since pumping started, 

» 

t = time in days since pumping stopped 
and 

3 

Q = rate of recharge = rate of discharge, in m /day. 

’ 2 
Uhen S and S are constant and equal and u = (r S/4Tt) is 

sufficiently small, equation 5.14 can be written as 


i s 


2.30 Q t 

log — 


(5.15) 


4 Tt T t 

For the data obtained in observation well or in the 

> f 

abstraction well, s vs t/t is plotted on semilog paper 

f 

(with t/t on logarithmic scale) and a straight line is 
fitted through the plotted points. The slope of the straight 
line is given by 2.30Q/4nT. The transmissivity (T) is 
obtained from the expression 


2.30 Q 

T = 

4 TT ^ S 

where is the drawdown for unit log cycle (of 

» ’ 

measured on s versus t/t plots. 


.. (5.16) 

» 

t/t ) as 
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5 . 3.6 


for 


Storativi ty 


Spsciflc Capacity Equation 
Determination 

Ualton (1970) reported an expression for specific 
capacity as related to transmissivity (T) and storativity (S) 
for a well pumped at a constant rate in a homogenous 
isotropic non-leaky artesian aquifer infinite in aerial 
extent. According to him, the specific capacity can be 
written as 


Q 


T 


264 


Tt 


log 


2693 r S 
w 


65.5 


where 

Q/s = specific capacity, in gpm/ft, 

Q = discharge in gpffl, 

s = drawdown, in feet, 

T = transmissivity in gpd/ft,. 

S = storativity, fraction, 

r^ = normal radius of well in feet, 

and 

t = time after pumping started, in minutes. 

If the following field units are taken, 

3 

Q/s in m /day/m 

Q in m^/day 

s in m 

2 

T in m /day 
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(5.19) 



in m 


w 


t in days 

then the equation 5.19 can be written as 
Q 26.3 X 10‘^ X T 


T t 

log - 7.71 


. (5.20) 


r-’w S 


Rearranging the terms, 


1 

log 

S 


26.3 X 10 


X T X £ 


Q 


+ 7.71 


log 


T t 


r 


* 

w 


(5.21) 

Thus the storativity of an aquifer can be obtained from 
specific capacity data. This procedure was adopted in the 
present study since no observation well facility was 
available. For this purpose, the values of Transmissivity (T) 
as obtained with Jacob’s (1950) method have been used. 
Although such estimates are not very accurate these were 
useful in characterizing the nature of the aquifer. 

5.4 PUMP TEST DATA ANALYSIS 

Pumptests have been conducted on selected tube 

wells located in the productive aquifer zone delineated 
within the three limestone members. The transmissivity and 
storativity values were estimated on the basis of time- 

drawdown and time-recovery data. 




Pmnptest results on j^lls located within the delineated zone 
5.4.1 Nipania ineinber 
1. Bhatapara 

The well is situated near Bhatapara Railway station. 
This well with a diameter of 0.2 m penetrates to a depth of 
6 7.07 ro . , through 42.60 m of aqui fer horizon. 

Q = 468.78 m^ / day 

t = 0.208 days 

Depth = 67.073 m 
^ s = 0 . 853 m 
^s’ = 0.749 m 

(a) Jacob’s Method using drawdown data, [Fig 5.4] 

2.3 X Q 

T = 

4 IT ^ s 

2.3 X 468.78 

= = 100.63 m*/<3ay 

4 X TT X 0.853 


(b) Recovery test [Fig. 5.5] 

2.3 X Q 


4 X IT Xi. s' 

2.3 X 468.78 
4 X TT X 0.749 


■ av 


= 107.62 m* /day 


(c) Storativity 


-7 


= 114.61 m’/ day 


26.3 X 10 X T X s 


Log 


+ 7.71 


log 


Tt 


w 


S 


Q 






\n . 

/ 



26.3X10 ^X107 .62X1.75 107.62X0.208 

+ 7.71 - log 


468.78 


(0.203)* 


log = 4.975 

S 


,-5 


S = 1.059 X 10 
2. Tilda 

Q = 645.04 /day 
t = 0.229 days 
Depth = 60 m 

= 0.20 in 
^ s = 0.82 in 
^ s ’ = 0.80 in 

(a) Jacob’s Method using drawdown data, [Fig. 5. 6 ] 

2.3 X Q 


T = 


4 T 14 s 

2.3 X 545.04 
4 X n X 0.82 

(b) Recovery test [Fig. 5.7] 

2.3 X Q 


121.71 in*/day 


T = 


4X n X A s’ 


2.3 X 545.04 
4 X TT X 0.80 
Tav = 12 3.20 in* /day 


= 124.69 in*/ day 
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FIG. 5-6 TIME- DRAWDOWN GRAPH FOR WELL AT TILDA 






AT TILDA 




(c) Storativity 


1 26.3 X 10 X T X s Tt 

Log = — + 7.71 - log ^ 

s o rJ 


26.3X10"^X123.2X1.64fi 


545.04 


+ 7.71 


log 


123.2X0.229 

(0.203)* 


1 

log = 4.874 

S 

S =1.33X10 


3 . Murpar 

The well is located in the vicinity of the post office. 
Q = 408.48 /day 
t = 0 . 2569 days 
Depth = 60 in 

= 0 . 1 5 m 
^s = 0.725in 
^ s ’ = 0 . 7 0 

(a) Jacob’s method using drawdown data, [Fig. 5. 8] 

2.3 X Q 

T = 

4 TIZ_ S 

2.3 X 408.78 

= 103.25 m’ /day 

4 X Ti X 0.725 

(b) Recovery test [Fig. 5. 9] 

2.3 X Q 

T = 

4X-n 
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^AWDOWN GRAPH FOR WELL AT MURPAR 








2.3 X 408.78 
4 X TT X 0 . 70 
Tg^-y = 105.09 m^/day 
(c) Storativity 


= 106.93 m /day 


Log 


26.3 X 10 ^ X T X s 


Q 


Tt 


+ 7.71 - log 


■ w 


26 . 3X10~^X105.09X1 .75 105.09X0.2569 

= + 7.71 - log 

408.78 (0.15)* 

1 

log = 3.0798 

S 

S = 8.33 X 10“^ 

5.4.2 Pendri Member 

In this member four pumping sites were taken in for analysis. 
1. Bitkuli 

3 . , 

Q = 590.46 ro /day 
t = 0.229 days 
ry = 0.15 
Depth = 60.97 m 
a = 0.53 m 
^ s = 0 . 52 5 m 

length of the slotted pipe = 7 m 

(a) Jacob’s Method using drawdown data, [Fig. 5. 10] 

2.3 X Q 

T ' = — 

A MAS 
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AT BlTKULl 



2.3 X 590.46 


4 X Ti X 0.53 

Cb) Recovery t est [ Fig . 5 . 1 1 ] 


203.90 m* /day 


2.3 X Q 

T = 

4 X -n Xas’ 

2.3 X 590.46 

4 X TT X 0.525 

Tav = 204.874 m* /day 


= 205.849 in’/ day 


Cc) Storativity 

1 26.3 X lO"^ XT X s 

Log — = + 7.71 - log 

S Q 


Tt 


2 


26 . 3X1 0~^X204 .874X1 .25 
590.46 

1 

log = 4.39 

S 

S = 4.05 X lO'^ 

2 . Balnsa 

Q = 13626 in^ /day 
t = 0.25 days 
Depth = 60.9 m 
Ty = 0 . 15 m 
A s = 2 . 2 in 
As’=2.12«i 


204.874X0.229 

7.71 - log 

CO. 15)* 
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AT BAINSA 




(a) Jacob’s Method using drawdown data, [Fig. 5. 12} 


2.3 X Q 

T = 

4 Tr _s s 

2.3 X 13(526 

4 X V X 2.2 

(b) Recovery test [Fig. 5. 13] 


1134 m* /day 


2.3 X Q 

T = 

4 X Tf X As’ 

2.3 X 13626 

= = 1176.98 m‘/ day 

4 X IT X 2.12 

Tav “ 1155.49 in* /day 

(c) Storativity 

1 26.3 X lO"^ XT X s 

Log — = + 7.71 - log 

S Q 


Tt 


• w 


2 


26.3X10 ^X1155.49X4.42 1155.49X0.25 

= + 7.71 - log 

13626 (0.15)* 

1 

log = 3.6 0 

S 

S = 2.50 X lO'^ 

3. Bhatbhera 

3 

Q = 610.44 m /day 
t = 0.26 days 
Depth = 59 . 45 m 

Ty = 0 . 20 m 

s = 0 . 65 m 
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0.65 m 


\ 


s’ 


(a) Jacob’s Method using drawdown data, [Fig. 5. 14] 

2.3 X Q 

T = 

4 TT ^ S 


2 .3 X 610.44 

= 171.97 in*/day 

4 X Ti X 0.65 


(b) Recovery test [Fig. 5. 15] 


2.3 X Q 

T = 

4 X TT X A.S’ 

2.3 X 610.44 

= = 171.97 m*/ day 

4 X TT X 0.65 


Tav = 171.97 m* /day 

(c) Storativity 


Log 


1 

S 


26.3 X 10~^ XT X s 

Q 


+ 7.71 


log 



26.3X10 ^X171 . 97X1 . 53 171.97X0.26 

= + 7.71 - log 

610.44 (0.203)* 

1 

log '■■■ =4.67 

S 

S = 2.1 X lO"^ 




4. Arjuni 

Q = 772.14 /day 
t = 0.236 days 
Depth = 60.97 m 


Ttj = 0.15 in 


^s = 0.975 m 
^ s ’ = 0.92 in 

(a) Jacob’s Method using drawdown data, [Fig. 5. 16] 


2 . 3 X Q 


4 TI s 


2.3 X 772.14 


4 X -n X 0.975 


144.99 in* /day 


Cb) Recovery test, [Fig. 5. 17] 

2.3 X Q 


4 TI s’ 


2.3 X 772.14 


4 X IT X 0.92 


Tov = 149.34 m /day 


(c) Storativity 


= 153.69 in* /day 


26 .3 X 10 X T X s 


Log 


+ 7.71 - log 


26. 3X10" X149.34X1.75 


772.14 


+7.71 - log 


149.34X0.236 


(0.152)* 


= 4.515 
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ECOVERY GRAPH FOR A WELL AT ARJUNl 


S = 3.05 X 10 

5.4.3 Newari Member 
1 . Khaira 

Q = 6540.48 in^ /day 
t = 0.104 days 
Depth = 56.40 m 
r ” 0.15 lu 
^ s = 1 . 7 m 
^ s ’ = 1.67 m 

(a) Jacob’s Method using drawdown data, [Fig. 5. 18] 

2.3 X Q 

T = 

4 TT A s 


2.3 X 6540.48 

= = 704.17 m’/day 

4 X TT X 1.70 

(b) Recovery test, [Fig. 5. 19] 

2.3 X Q 

T = 

4 TT / S ’ 


2.3 X 6540.48 


4 X TT X 1.67 


717.1 m'/day 


■av 


= 710.63 m /day 


(c) Storativity 


26.3 X 10 ^ X T X s 


Tt 


Log 


4 7.71 - log 


■ w 
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AT KHAIRA 



+ 7.71 


710.63X0.104 


26.3X10 ^X710 . 63X1 . 75 
6540.48 


- lOfi 1- 

(0.15)* 


1 

lofi = 4.193 

S 

S = 6.404 X lO”^ 

5.4.4 Performance Characteristics of the Wells Located 
Outside the Delineated Zone 

To enable comparison of the performance of the 
formations as aquifers, tests have also been conducted on 
five wells located outside the identified zones. The 
transmissivity estimates are as follows : 


Transmissivity (T) m* /day 

Member 

Well location Time-drawdown 

Time-r ecovery 

Average 

NIPANIA 

KUKRACHUDA 

51.3 

56.3 

53.8 

NIPANIA 

DHURRABANDH 

58.8 

64.6 

61.7 

PENDRI 

PUR ENA 

65.0 

69.0 

67.0 

NEUARI 

BHATGAON 

141.0 

149.0 

145.0 


As can we seen from the results, the Pendri member has 

2 

re£istered high transmissivity values (of around 1100 m /day) 
in the Bainsa area. It may be indicated that the digital 
processing and subsequent ground checks have revealed that 

the limestone member in the area to be closely fractured and 
highly cavernous (Fig.^'M. The trahsmissivity estimates in 
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e o. er ree localities within this member range between 

144 and 204 m /day. Present investigat ions revealed 

prominence of fractures in the localities. In contrast, the 

transmissivity value for Pendri member registered at Purena 
2 

is around 61 m /day. This locality is out side the delineated 
zone and has relatively compact limestone devoid of cavities 
and fractures. 

For the three wells in the Nipania formation within the 

delineated zone the transmissivity estimates are in the same 

2 

range of around 100 to 110 m /day. However, outside this 

zone, the formation has registered comparatively lower values 

2 

of transmissivity of around 53 m /day in the Kukrachunda area 

2 

and around 61 m /day in the Dhurrabandh area. 

In the Newari formation, while the transmissivity 

2 

estimate was around 704 m /day in the Khaira area, the value 

2 

drastically is low (of around 145 m /day) around Bhatgaon due 
to the compact natures of the limestone. The storativity 
values in all the cases are in the range for confined 
aqui f ers . 
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CHAPTER 6 


SUMMARY AND CONCLUSIONS 

In the carbonate terrains, delineation of aquifer 
zones involves recognition and mapping in detail of the 
various features such as the fracture zones, sinkholes, 
caverns and collapsed zones, that aid in the ability of 
formation not only to store water but also to transmit the 
same. Lithology plays an important role in the selective 
dissolution of the formation resulting in solution features. 
Since many years, visual interpretation of the imageries in 
remote sensing has been carried out in delineating the 
morphological and structural features. Thermal infrared 

imageries have also been used in the visual interpretation 
However, digital image processing enables an accurate 

interpretation as the digital data is involved. A combination 
of image processing techniques using MSS data has facilitated 
a clear demarcation of the structural and morphological 
features in the karstified zones. Uhile band ratioing gives 
the effects of slope change within the valley region, 
stretching improves the overall resolution to enable a better 
interpretation. Edge enhancement with proper filter provides 
not only the enhancement but also creates an edge for all the 
subdued caverns and collapsed zones in a carbonate terrain. 
Selective stretching results in the enhancement of a 

particular area of interest, in the' present case caverns have 

% 



been enhanced through this technique. All these techniques 
have been utilized in locating aquifer zones vithin the 
limestone terrain. 

Central, part of the Chhatt is,garh basin is taken up 
as a case study for the application of methodology for 
digital image processing developed by the author during the 
course of the present work. Located in Hadhya Pradesh, the 
study area covers parts of Raipur, Baloda Bazar and Bilaspur 
districts. The formations in the area belong to the Raipur 
series. This series essentially, is a repetitive sequence of 
limestone and shale with 5 cycles occurring in the basin. A 
sub-tidal to tidal environment has been inferred for the 
deposition of these sediments (Schnitzer , 1974) . The formation 
of the shale has been attributed to release of fine mud 
through minor uplift of landmass by epirogenic forces. 
Predominance of limestone in the sequence is indicative of 
prevalence of warm temperatures and chemical weathering in 
the provenance. Variations in the lithology of shales 
revealed changes in the oxidizing and reducing conditions in 
environment during deposition. The formations within the 
study area pertain to the third and fourth cycles. The 
Gunderdehi shale (of the third cycle) is overlain by the 
limestone members (Nevari, Pendri and Nipania) and the 
Tarenga shale belonging to the 4th cycle. The Pendri and 
Newari limestones are .stromatolitic:- Variations in lithology 



recogoi t iO'jn of 


in tihose limestone members have facilitated 
the same on the images. Hitherto, no systematic 
classification of this carbonate terrain has been attempted 
by the earlier workers in terms of morphological feattires, 
drainage characteristics and the distribut ion of the aquifer 
zones. 

In the present work the data on tape from Landsat-4 

was windowed on OMEGA system and this data has been converted 

into Dec mode for image classi f icat ion . Initially the window 
data when viewed on PDF 11/24 system was found to possess a 
skew effect as consequence of which the generated picture of 
the area did not match with the toposheet of the same area. 

This skewness effect, due to Earth's rotation has been 

corrected using nearest neighbourhood method. Software 
necessary for all corrections has been developed in present 
work. The corrected window data has been processed through 
several steps of image processing in sequence. Band ratioing 
was adopted as a first step and it has been found that, of 
all band combinations , the ratio of bands 2 and 4 yields 
better results in enhancing the slope of the river valleys. 
In the paeudocoloured technique, the corrected window of the 
band data has been pseudocoloured by intensifying the hue of 
the colour, which gives better results of surfacial feature 
identification. Tonal contrast is obtained for variations of 
reflectance all over the terrain. The sinkholes appeared hazy 
in the ratio pictures. Stretching has been done, on the ratio 
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picture data (based on the distribution of the reflectance 
values) to obtain better resolution of the features. The 
technique involves plotting histogram for the data initially 
covering range 1.0 - 2.5 (reflectance values) which has been 
later stretched to the entire dynamic medium. This has 
resulted in the improvement of the intensity and enhancement 
of the sinkholes which were otherwise observed to be hazy in 
the initial ratio pictures. Further enhancement of caverns 
was achieved by subjecting the data to selective stretching 
with the help of the reflectance value ranges of the 
cavernous areas. The study area has been classified on the 
basis of tonal characteristics into three major morphological 
units- the upland tract, the valley tract and the alluvial 
tract- The upland tract is characterized by its flat 
topography. This tract has been further subdivided into two 
units, based on the variation in reflectance values. The 
valley tract is charact erized by typical pattern of gully 
erosion as a consequence of the dissolution of limestone. The 
alluvial tract, situated in the northeast corner of the area, 
is characterized by a gray tonal contrast. 

The edge enhancement technique used in the present 
study has facilitated the delineation of the edges of 
collapsed and sinkhole-prone zones. For this, the noise has 
been smoothened through a correction using 3 by 3 filter 
initially. However, since this procedure has not yielded all 


i63 




MAP INDICATING POTENTIAL AQUIFER ZONES DELINEATED 
THROUGH IMAGE PROCESSING IN THE STUDY AREA 






the details. 


^ 5 by 5 filter was developed by the author 
which has provided sharp boundaries wherever the reflectance 
values change suddenly. This technique has enabled distinct 
demarcation of the collapsed zones in the study area. 
Subsequently the same data has been subjected to convolution 
operation with gradient filter. This procedure has helped in 
highlighting all the lineaments. Detailed mathematical 
treatment of the procedures has been presented together with 
the relevant soft ware developed by the author. The 5 by 5 
filter has also been used in enhancing the dense vegetational 
areas within the study region through histogram equalization 
of the window data. 

Image classification is done by taking the a priori 
probability of features into consideration. The data has been 
checked with a numerical separation technique using the 
training data set and a line printer map of the study area 
was prepared. 

Uithin the limestone terrain, productive aquifer 
zones [Fig. 6.1] have been delineated in the study area using 
the image processing techniques. Pumptests have been 
conducted on a few wells located in the three limestone 
members within this zone and the estimates of transmissivity 
and storativity were obtained on the basis of time drawdown 
and time-recovery data. The very low values of storativity in 
the range of 10~^ to 10~^ reflect the confined nature of the 
aquifer. The Fendri limestone has registered high 
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transmissivity values ranginfi between 200 and 1155 m^/day 
within the aquifer zone depending on the presence of the 
fractures, sinkholes and the cavernous features. In the 
localities where the fractures are present but caverns are 
not seen a.s in Nipania and Nevari inerobers , the 

t ransmiss i vi t y is around 124 to 700 m^/day. Puinptests carried 
out on veils located outside the delineated aquifer zone 
within these three limestone members have revealed relatively 
low transmissivity values of the order 53 to 61 m^/day. Thus 
the productive aquifer zone, as delineated on the basis of 
digital image processing in the present work has been proved 
for ground water development, through aquifer evaluation 
t ests . 

Uater level contour maps for pre-monsoon and post- 
monsoon periods have enabled the identification of recharge 
areas. The general flow pattern is towards the center of the 
study area where the caverns and sinkholes in the Pendri 
member act as storage spaces for groundwater. The quality of 
groundwater has been observed to vary in flow direction. The 
TDS content and the hardness have registered an increase in 


the flow direction. 
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APPENDIX I 


EVALUATION OF REHOTE SENSORS AS HYDROLOGIC TOOLS 

(ROBINOV 1968) 


General Comments on Potential value 
and use of Remote-Sensor System for 
Remote-Sensor System Hydrologic Studies 


Panchromatic Panchromatic photography is the most 
Photography widely used remote sensing technique 

because of its availability and rela- 
tively low cost. Interpretive techni- 
ques are well developed and formal 
training of its use is available. 


Mult ispectral 
Photography 


Multispectral photography interpreta- 
tion requires a background of spectral 
signature studies of terrain and water 
features that have not yet been made. 
Data returns from multispectral 
systems may be so voluminous that they 
can not be readily interpreted. 
Little work has been done on 
interpretation for hydrologic 


purposes . 


Infrared 

Photography 


Colour Photography 


Colour Infrared 


Infrared photography is primarily of 
value in mapping drainage features and 
shorelines- The w^ater is always black 
in positive print. Some vege-tation 
characteristics are discernible. Its 
most valuable use is as an adjunct to, 
but not a replacement for, standard 
aerial photography - 

Colour photography, in spite ol its 
built-in spectral redundancy, promises 
to be a major tool of hydrologist in 
many special fields and is suffic- 
iently better for recognition of 
significant hydrologic features that 
it may replace panchromatic photo- 


graphy 

for many uses. 

The interpr- 

station 

capabi 1 i ty 

of 

the potential 

operational hydrologic 

users of 

colour 

photography 

must be greatly 

increased. Hethods 

for 

spectral and 

density 

extraction 

of data are being 


developed . 

Colour-infrared photography may , be 
superior to standard phtography in 
some respects.' It shows di f f erences .in 



vegetation more clearly and , provides' a 


slightly higher contrast on water 
surfaces. Its general superiority to 
standard colour photography has yet to 
be proved but it may be highly useful 
and its worthy of much additional 
research . 

Infrared Radiometry Infrared radiometry is very useful for 

sequential measur ements of changes in 
land and water surface temperatures 
because it is a simple measurement 
technique and data reduction is 
simpler than for infrared imagery. 
Radiometry is routinely used for 
periodic surveys of near shore oceanic 
areas. 

Infrared Imagery Infrared imagery has shown its value 

as a tool for measuring water surface 
temperature and as a means of 
qualitatively dif f erentiating some 
terrestrial features. The lack of a 
simple means of determining eroissivity 
hampers its quantitative usefulness. 
Analytical techniques for proper use 
of the reduced data need to be 


developed . 



Radar Imagery 


Side-looking airborne radar has an all 


weather capability for coverage of 
large areas. Its ability to penetrate 
foliage and accentuate topographic 
features enhances its value, Uater 
surfaces are excellent reflectors of 
Hiicrovaves, resulting in a uniform 
black-tone image. For these reasons 
stream drainage systems and water 
surfaces are easy to identify. The 
black -tone precludes measuring the 
physical, chemical or biologic 
characteristics of water. Radar may be 
of value in terrain analysis for 
ground water exploration. 

Hicrovave Radiometry Passive-microwave sensors measure the' 
and Imagery brightness temperature of terrain and 

water surfaces. Spatial resolution is 
lower than infrared systems but 
radiance is directly proportional to 
temperature. Probably will find 
greatest application in oceanic and 


snow-f ield mapping. 



APPENDIX II 


FEATURES THAT ARE IMPORTANT FOR MAPPING CONSOLIDATED ROCK 
AQUIFERS ON LANDSAT MSS IMAGES (FROM MOORE, 1982) 

ROCK TYPE : 

1. Landform; topographic relief; erosional characteristics. 

2. Outcrop patt ern--banded pattern for sedimentary rocks ;lobate 
outline for basalt flows; arcuate or faulted boundaries 
for igneous intrusions. 

3. Drainage pattern and drainage texture (density). 

4. Fracture type and symmetry (as indicated by lineaments). 

5. Topographic position; steepness of slope; presence or 
absence of erosional terraces on hillsides. 

6. Shape of ridgeline and plateau edge. 

7. Type and density of vegetation cover; landuse. 

8. Presence of sinkholes (visible only where occupied by water 
or distinctive vegetation). 

9. Tone or hue; image texture. 

FOLDS : 

1. Ridge and valley topography; mountain, dome, and elliptical 
hill; short, parallel and intersecting ridge crests with V- 
shaped sections (complexly folded metamorphic rocks); cuesta 
or hogback. 

2. Flatirons on dip slope and irregular topography on back 



slope; parallel to dendritic drainage on dip slope and 
angulate or trellis drainage on back slope; uniform 
distribution of vegetation on dip slope and vegetation 
banding parallel to ridge crests on back slope. 

3. Banded out crop pattern not related to topography. U-shaped 
to V-shaped map pattern of these ridges. 

4. Trellis, radial, annular, or centripetal drainage pattern; 
partly developed pattern of these types. 

5. Major deflection in stream channel; change in meander wave- 
length or change from meandering to straight or braided 

channel . 

6. Assymmetric drainage; channel not centered between drainage 
divides. 

LINEAMENTS : 

1. Continuous and linear stream channel or valley; 

discontinuous but straight and aligned valleys, draws, 

# 

swags, and gaps. 

2. Elongate or aligned lakes, large sinkholes, and valvanoes. 

3. Identical or opposite deflections (such as doglegs) in 

adjacent stream channels or valleys; alignment of nearby 
tributaries and tributary junctions. 

4. Elongate or aligned patterns of natural vegetation; thin 

strip of relatively open (may be right-of-way) or dense 

vegetation. 

5. Alignment of dark or light soil tones. 
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time in minutes that must elapse for straight-line 

METHOD TO BE APPLICABLE (After S.H. SIDDIQUI ) 
time in MINUTES FOR r. , = 0 . 07 ^ tn 




CHEMICAL ANLYSIS-RESULTS OF WATER SAMPLES 
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******************************************************* 

* This Program converts the latitudes and longitudes * 

* of a given point on the Earth’s surface to its line * 

* number and pixel number on the Imagery after * 

* conversion to conical orthomorphic coordinates. The * 

* constants in this program are valid only for the * 

* Chhattisgarh area. * 

******************************************************* 

IMPLICIT DOUBLE PRECISION(A-H , K-Z ) 

DIMENSION LAT(150) ,LONG(150) .LATRAD(150) 

LONGO=81 . 0 

LAT0=21 .0*3.1415926/180.0 
A=6377277 . 6 
ELLIP=1 .0/300.8 
ELLIP=ELLIP*ELLIP 
SIN=DSINCLA-T0) 

SINS0=SIN*SIN 
DEN0M=1 .0-ELLIP*SINSQ 
N0=A/DSQRT(DEN0M) 

R0=N0*CC1 .0-ELLIP)/DENOM) 

P=N0*(DCOSCLAT0)/DSINCLAT0)) 

RM=DS0RT(R0*N0) 

URITEC22,18) 

READ(21,*)II 
DO 100 1=1,11 

READC21 , * )DEG1 ,MIN1 ,SEC1 ,DEG2 ,MIN2 .SEC2 

LAT( I)=DEG2+MIN2 /60 . 0+SEC2/3600 . 0 

L0NG( I)=DEG1+MIN1/60.0+SEC1/3600.0 

LATRADC I )=LATC I) *3. 1415926/180.0 

L0NGDF=LONGC I )-LONG0 

GAMMA=LONGDF*SIN 

GAMMA=GAMMA*3. 1415926/180.0 

M=RM* (LATRADC I) -LATOl 

M1=M**3/(6.0*RO*N01 

TAN=DSINCLAT0)/DCOSCLAT0) 

M2=M1*(M*TAN/(4.0*N0)) 

M11=M+M1+M2 

X=(P-Mll)*DSIN(GAMMA) 



GAHMA=GAHnA/2 .0 

Y=I11 1 + X * ( DS I N ( G AHHA ) / DCOS ( GAMHA ) ) 

Al= -15.37834371515483 

Bl= 54.28787394037448 

C1=X - 75235.19856468876090 

A2= -76.71755874202162 

B2= -10.81379684548759 

C2=Y - 184572.45036062989400 

X1 = CB2*C1-B1*C2 )/ (A1*B2-A2*B1) 

ILIN=Xl+0.5 

ILIN1=ILIN*4 

Y2 = CA2*C1-A1*C2)/(A2*B1-A1*B2) 

IPIX=Y2+0.5 

URITEC22 , 1001 )LAT(I ) ,L0NG(I ) , ILIN, IPIX , I LINl 

FORMAT (7X, 60 ClH-)//9X, 'LATITUDE' ,6X, 'LONGITUDE' ,4X, ’LINE NO’ , 
13X, 'PIXEL NO’ ,3X, 'RECORD NO * , / /7X . 60 ( IH- ) / ) 
FORMATC6X,F12.8,3X,F12.8.4X, I4,6X,I4.7X,I4)' 

CONTINUE 

F0RMAT(5X, ’LATITUDE=’ ,F23.16,5X, ’ LONGITUDES= ' ,F23.16/) 
FORMAT(5X, 79(1H=)//3X, 'LATITUDE' ,5X, 'LONGITUDE' , / /5X , 7 9 ( IH- ) / 'I 
STOP 
END 
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* 

* 

* 

* 

* 

* 

* 


* 

INTERACTIVE PROGRAM FOR HISTOGRAM * 

Reads a full-resolution (256 in) histogram from the * 
selected image refresh memory, then plots the results* 
in the selected overlay plane, enabling overlay * 
visibility. The process may be repeated indefinitely.* 

* 


* BYiC.D.MURTY Date:24-5-85 * 


EXTERNAL lAOVLY, ICSR, IBPMA, lOPOP 

DIMENSION HIST(256), IHIST(256), ICHARCIOOO) 

EQUIVALENCE (HIST(l), ICHAR(l)), (IHISTCl), ICHAR(513)) 

JAOVLY = lADDR (lAOVLY) 

JCR = lADDR (ICSR) 

JBPMA = lADDR (IBPMA) 

JOPOP = lADDR (lOPOP) 

CALL IP5INI 


CALL SETUP 

CALL I POKE (JOPOP, ”41503) 

DO 10 I = 1, 1000 
ICHAR(I) = 0 
CONTINUE 

ENCODE (7, 900, ICHAR(l)) 

FORMAT ( ’ CHANNEL ’ ) 

ENCODE 7, 910, ICHAR(121)) 

FORMAT ('MAXIMUM') 

CALL BLKMOV (ICHAR, lAOVLY, 1000) 

CALL IPOKE (JCSR, lOR (IPEEK (JCSR), ”20000)) 
TYPE 1000 

FORMAT (’$Channel to be histogrammed (0-4): ’) 

ACCEPT 1010, NCHAN 
FORMAT (II) 

IF (NCHAN -GT. 4) GO TO 20 
MHS = 2**NCHAN 
TYPE 1020 

FORMAT (’$Channel for plot, bit plane: ) 

ACCEPT 1030, IPLCH, IBIT 
FORMAT (212) 



40 


50 


60 


920 


-EO. 

1) 

NBIT 

.LE. 

1) 

CALL 

.GE. 

2) 

CALL 

HIST) 




= ISUAB (MBIT) 

IPOKE (JBPHA, NBIT) 

I POKE (JBPHA 2, NBIT) 


IF (IPLCH -LT. 0 . OR . I P LCH . GT . 4 ) GO TO 30 

IF (IBIT .LT. 0 .OR. IBIT .GT. 7) GO TO V 

IF (IPLCH .EQ. 4) CALL OVLINI 

MPLCH = 2**IPLCH 

mbit = 2**IBIT 

NBIT = MBIT 

CALL OUTNBL (MPLCH) 

IF (MOD (IPLCH, 2) 

IF (MOD (IPLCH, 4) 

IF (MOD (IPLCH, 4) 

CALL FMOPER 
CALL GTHSTA (MHS, 

RMAX = 0.0 
DO 40 I = 1, 256 

IF (HI ST (I) .GT. RMAX) RMAX = HIST(I) 
CONTINUE 

SCALE = 5 11. 9 9 /RMAX 
DO 50 I = 1, 256 

IHIST(I) = SCALE*HIST(I) 

CONTINUE 
DO 60 I =, 1 , 256 

IX = (I - 1)*2 
CALL VECTOR (IPLCH, IX, 0 
CONTINUE 

CALL IPOKE (JCSR, lOR (IPEEK (JCSR) 

ENODE (1, 1010, ICHAR(769)) NCHAN 
ENCODE (8, 920, ICHAR(770)) RMAX 
FORMAT (1X,F7.0) 

CALL IPOKEB (JAOVLY + 83, ICHAR(769)) 

CALL BLKMVA (ICHAR(770), JAOVLY + 320, 4) 

GO TO 20 
END 


IX. IHIST(I), MBIT, MBIT) 


‘40)) 


******************************************************* 

* SUBROUTINE GTHSTA (MEMORY, HIST) 

******************************************************* 


* 

* 
t * 


DIMNNSION HIST(256), IHIST(512) 

MEMSEL = MEMORY 

CALL RDHSTA (MEMSEL, IHIST) 

DO 10 I = 1, 256 
K = 1*2 

j = K - 1 

IF (IHIST(J) -LT. 0) GO TO 5 
HIST(I) = 65536. 0*IHIST(K) + IHIST(J) 

GO TO 10 

HIST(I) - 65536.0*(IEIST(K) +1) + IHIST(J) 


5 



10 


CONTINUE 

RETURN 

END 



c - * 

c INTERACTIVE PROGRAM FOR BAND RATIOING * 

C (Using some System routines) * 

C (1) the memory channel containing the dividend, * 

c (2) the memory channel containing the divisor and * 

C (3) the memory channel in which the remainder will * 

C appear. * 

C The quotient overwrites the the dividend. * 

C * 

C BY : C.D.MURTY DATE:24.5.a6 * 

Q* *1c ** -kit t **** -k ****** -k -k **** -k **** -k 1c *********** -k *************** * 

c 

EXTERNAL ITSTOP, lOPOP, lODPSH 

CALL IP5INI ! Set up mapping to IP-5000 

JTSTOP = lADDR (ITSTOP) 

JOPOP = lADDR (lOPOP) 

JODPSH = lADDR (lODPSH) 

1 TYPE 1010 

1010 FORMAT (’JChannels for dividend, divisor, 

remainder : ’ ) 

ACCEPT 1020, MSQ, IDV, IRM 
1020 FORMAT (312) 

IF (MSQ .LT. 0 .OR. MSQ .GT. 3) GO TO 1 

IF (IDV .LT. 0 .OR. IDV .GT. 3) GO TO 1 

IF (IRM .LT. 0 -OR. IRM .GT. 3) GO TO 1 

IF (MSO .EQ. IDV .OR. MSQ -EQ. IRK) GO TO 1 

IF (IDV .EQ. IRM) GO TO 1 

MSQ2 = 2**MSQ 

IDV2 = 2** IDV 

IRM2 = 2**IRM 

CALL SETUP 

CALL INSEL (MSQ2, 2, 0) 

CALL IPOKE (JTSTOP, ”20) 

CALL IPOKE (JOPOP, ”41503) 

CALL IPOKE (JODPSH, IRM2*”400 + ”331) 

CALL OUTNBL (IRM2 + MS02 ) 

CALL FMOPER 

CALL INSEL (IRM2, 10, MS02, 4, IDV2, 1. 34) 

CALL IPOKE (JTSTOP, ”2006) 

CALL IPOKE (JOPOP. ”50106) 

CALL IPOKE (JODPSH, IRM2*”400 + ”331) 

CALL OUTNBL (MS02 + IRM2) 

DO 10 ISTEP =1,7 
CALL FMOPER 
1 CONTINUE 

CALL IPOKE (JODPSH, IRM2*”400) 



CALL FMOPER 

STOP 

END 



00100 c******************************************************** 


00200 C c 
00300 C THIS PROGRAM USES 5X5 CONVOLUTION OPERATOR FOR C 
00400 C FINDING EDGES ON IMAGE DATA STORED IN SPECIFIED C 
00500 C CHANNEL. IT CALLS SUBROUTINE ’CONVOL' FOR PERFORMING C 
00600 C CONVOLUTION. C 
-00700 C C 
00800 C AUTHORzC.D.MURTY DATE : 15- JAN-1987 C 


00900 C*************^****’^******^*********’‘^****5*r*3*f5fct#*****il:**** 

01000 
01100 
01200 
01300 
01400 

01500 
01600 
01700 

01800 1030 

01900 
02000 

02100 1040 

02200 
02300 
02400 
02500 


02600 

c— -- 



02700 

02800 

c 

SUBROUTINE CONVOL 

C 

c 

SUBROUTINE CONVOL C FACTOR , ISOUR2 , ITARG2) 


02900 

03000 


DIMENSION FACTOR(5,5), ITTC256) 


03100 


IMPLICIT INTEGER CA-2) 


03200 


LOCAL IPOS 


03300 


EXTERNAL I POP , lODPSH , I TSTOP 


03400 


DATA ITT/256*0/ 


03500 


JOPOP=IADDR(IOPOP) 


03600 


JODPSH=IADDR( lODPSH) 


03700 


JTSTOP=IADDR(ITSTOP) 


03800 


CALL SETUP 


03900 


CALL 2EROM(ITARG2) 


04000 


CALL NBLITT(ISOUR2, 2) 


04100 


CALL INSEL(ISOUR2,4,ITARG2,8.0) 


04200 


CALL IPOKEC JODPSH, ”400*ITAG2+”300) 


04300 


CALL OUTNBLCITARG2) 


04400 


DO 80 1=1,5 



EXTERNAL I CSR , ICBCRG , IMAGE , ITSTOP , IPPOP , I ODPSH 
DIMENSION FACTOR (5, 5 I 
IMPLICIT INTEGER(A-2) 

DATA FACTOR /I, 0,2, 0,1, 0,4, 2, 4, 0,2. 2, -36, 2, 2, 0.4, 
2, 4, 0,1, 0,2, 0,1/ 

CALL IP5INI 
CALL SETUP 
TYPE 1030 

FORMAT ( ’ $SOURCE CHANNEL : ’ ) 

ACCEPT *,SCHAN 
TYPE 1040 

FORMAT ( ’ $TARGET CHANNEL : ’ ) 

ACCEPT *,TCHAN 

CALL CONVOL(FACTOR,2**SCHAN, 2**TCHAN) 

STOP 

END 



04500 


DO 70 J=l,5 

04600 


IF(FACTORCI , J) .EO.O.O)GOTO 70 

04700 


IXSHFT=I-5 

04800 


IYSHFT=J-5 

04900 


AFAC=IABSCFACTOR(I , J)> 

05000 


IPOS=.FALSE. 

05100 


IFCFACTORd , J) .GE.O.O)IPOS^.TRUE- 

05200 


DO 30 K=1,256 

05300 


ITT(K)=ISHFT(K-l,-5) 

05400 

30 

CONTINUE 

05500 


CALL READY 

05600 


CALL SCRLZMC IS0UR2 , IXSHFT, lYSHFT) 

05700 


CALL L0ADIT(IS0UR2, ITT) 

05800 


IF(IPOS)CALL IPOKECJOPOP, ”54531) 

05900 


IFC -NOT. IP0S)CALL IPOKECJOPOP , ”43106 ) 

06000 


DO 60 ITIME=1,AFAC 

06100 


CALL FMOPER 

06200 

6 0 

CONTINUE 

06300 

70 

CONTINUE 

06400 

80 

CONTINUE 

06500 


DO 90 1=1,128 

06600 


ITT(I)=(I-1) 

06700 

90 

CONTINUE 

06800 


DO 100 1=129,256 

07100 


ITT(I)=(257-I) 

07200 

100 

CONTINUE 

07300 


CALL READY 

07400 


CALL L0ADITCITARG2,ITT) 

07500 


CALL NBLITTC ITARG2 ,2) 

07600 


CALL IPOKECJOPOPP, ”47517) 

07700 


CALL FMOPER 

07800 


CALL SCRLZMC ISOUR2, 0,0) 

07900 


RETURN 

08000 


END 



